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1  FOREWORD 


In  this  project,  a  novel  control  technique,  referred  to  as  Closed  Loop  Vibrational  Control, 
has  been  developed  and  applied  to  the  problem  of  fuselage  vibrations  suppression  in  helicopter 
dynamics.  From  the  theoretical  standpoint,  the  technique  developed  is  applicable  to  systems 
where  the  control  input  enters  the  open  loop  dynamics  as  an  amplitude  of  a  periodic,  zero 
average  function,  and  this  amplitude  can  be  chosen  to  depend  on  the  system’s  outputs.  An 
example  of  such  a  system  is  the  helicopter  with  Higher  Harmonic  Control  (HHC)  where  periodic 
feathering  of  rotor  blades  around  a  fixed  pitch  angle  is  introduced  in  order  to  suppress  the 
fuselage  vibrations.  From  the  practical  standpoint,  the  technique  developed  is  useful  for  plants 
where  conflicting  control  objectives  must  be  achieved  with  an  insufficient  number  of  actuators. 
From  this  perspective,  the  technique  developed  is  b2ised  on  the  frequency  separation,  i.e.  the 
utilization  of  low  and  high  frequency  control  signals  so  that,  on  the  average,  all  control  objectives 
are  satisfied.  In  the  HHC  case,  this  frequency  separation  amounts  to  low  frequency  rotor  blades 
pitch  angle  control  to  ensure  the  desired  altitude  of  the  hovercraft  and  the  high  frequency  rotor 
blade  pitch  angle  control  to  suppress  the  fuselage  vibrations. 

For  systems  with  this  structure,  the  following  problems  have  been  solved  and  are  reported  in 
this  document: 

1.  Conditions  for  the  state  and  dynamic  output-feedback  stabilizability  by  closed  loop  vibra¬ 
tional  control  have  been  derived. 

2.  Pole  placement  capabilities  of  vibration2d  controllers  have  been  investigated. 

3.  Stability  robustness  of  closed  loop  vibrational  control  has  been  analyzed. 

4.  Youla-type  parametrization  of  closed  loop  vibrational  controllers  have  been  derived  and 
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utilized  for  the  design  purposes. 

5.  A  method  for  /f2“optinial  zeros  placement  has  been  developed. 

6.  The  results  obtained  have  been  applied  to  helicopter  vibration  suppression  problem  and  a 
technique  referred  to  as  Very  High  Harmonic  Control  (V'^HHC)  has  been  investigated. 

In  short,  the  main  result  can  be  formulated  as  follows;  A  novel  control  technique  has  been 
developed  and  its  utility  in  helicopter  vibrations  suppression  heis  been  demonstrated. 
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4  STATEMENT  OF  THE  PROBLEM 


The  goal  of  this  thesis  is  the  development  of  a  control  theory  for  a  class  of  dynamical  systems 


described  by  the  following  equations: 

i(0  = 

Ax{t)-i-  Bu{i)f{-)  . 

t 

(0.1) 

y{t)  = 

Cx(t) . 

fir)  = 

/(r  -h  T)  ,  r  7^  0  , 

~£n-r)<ir=0, 


0  <  c  <  1  , 

where  x  €  JR"  is  the  state,  y  €  iR  is  the  output,  u  £  JR  is  the  control,  f(t)  is  a  periodic, 
average  zero  scalar  function,  and  c  is  a  small  positive  parameter.  Stabilizability  properties  of 
system  (0.1)  with  state  and  output  feedback  are  analyzed  and  the  pole  placement  capabilities 
investigated.  A  characteristic  feature  of  system  (0.1)  is  that  the  control,  u,  enters  the  open 
loop  dynamics  as  an  amplitude  of  a  periodic,  zero  average  function.  Such  situations  arise  in 
a  number  of  applications  where  two  conflicting  control  goals  have  to  be  accomplished  by  a 
single  actuator.  For  instance,  in  the  helicopter  control  problem,  a  single  actuator  (the  blades’ 
pitch  angle)  is  used  to  ensure  both  the  desired  altitude  and  the  fuselage  vibration  suppression. 
These  goals  are  conflicting  in  the  sense  that  if  the  pitch  angle  is  chosen  to  ensure  the  desired 
2Jtitude,  the  fuselage  vibrations  are  not  suppressed;  if  the  pitch  angle  is  used  to  suppress  the 
fuselage  vibrations,  the  desired  altitude  is  not  attained.  In  order  to  accomplish  the  two  goeils 
simultaneously,  a  frequency  separation  approach  may  be  employed.  Specifically,  a  low  frequency 
control  may  be  used  to  stabilize  the  desired  altitude  and  a  high  frequency,  average  zero  control 
may  be  used  to  suppress  the  vibrations  without  compromising  the  first  goal.  When  the  control 
loop  is  closed  with  respect  to  the  low  frequency  control,  the  equations  have  the  form  of  system 
(0.1),  and  the  goal  is  to  choose  the  control,  u,  as  a  function  of  x  or  y  so  that  the  resulting  system 
has  the  desired  dynamical  properties. 

In  particular,  the  above  ideology  has  been  successfully  implemented  in  the  Higher  Harmonic 
Control  (HHC)  of  helicopters,  where  periodic  feathering  of  rotor  blades  around  a  fixed  pitch  angle 
is  introduced  in  order  to  suppress  the  fuselage  vibrations.  Helicopter  vibration  is  a  long  standing 
problem,  ^cent  experiments  [iMfi]  have  shown  that  HHC  may  lead  to  an  order  of  magnitude 
reduction  in  fuselage  vibrations.  The  primary  difficulty  in  implementation  of  the  HHC  systems  is 
the  complex  interaction  between  inertia,  structural  and  aerodynamical  hub  shears  and  moments 
for  HHC-equipped  helicopter  rotors.  These  interactions,  which  are  difficult  to  predict  due  to 
the  highly  complex  dynamics  of  most  helicopter  systems,  account  for  the  extreme  sensitivity  of 
HHC  efficacy  to  proper  magnitude  and  phasing  of  the  HHC  inputs.  Using  the  idea  of  closed 
loop  vibrational  control,  we  will  explore  alternative  means  of  suppressing  the  vibratory  airload 
without  the  phase  dependencies  in  the  input. 

It  is  well  known  that  unavoidable  discrepancies  between  mathematical  models  and  real-world 
systems  can  result  in  the  degradation  of  control  system  performance.  Thus,  wt  also  investigate 
the  property  of  stability  robustness  for  system  (0.1).  Both  synthesis  and  analysis  problems  are 
addressed.  In  the  synthesis  problem,  it  is  assumed  that  (0.1)  is  the  nominal  plant,  whereas  the 
true  plant  is  defined  by 

X  =  (A  +  AA)x  +  Buf  ,  (0.2) 
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where  A/l  is  the  perturbation  matrix.  Given  the  perturbation  matrix  A.4.  we  derive  a  condition 
which  guarantees  the  existence  of  a  controller 


e 


which  robustly  stabilizes  the  uncertain  system  (0.2).  In  the  analysis  problem,  we  determine  a 
bound  on  A/1  which  ensures  that  a  controller  (0.3)  which  stabilizes  the  nominal  plant  (0.1),  also 
stabilizes  the  perturbed  system  (0.2). 

We  also  investigate  the  parametrization  of  stabilizing  controllers  for  the  system  (0.1)  since 
such  a  characterization  could  be  quite  useful,  e.g.  to  satisfy  desired  performance  specifications 
[6].  We  first  consider  the  observer-based  output  controller,  Knomi  defined  by 

X  =  Ax  +  Buf  L{y  -  y)  ,  (0.4) 

K . 

u  =  — X  , 

y  =  Cx  , 

and  parametrize  all  stabilizing  controllers  for  (0.1)  in  the  class  of  rational  transfer  functions.  Next, 
the  parametrization  of  the  averaged  closed  loop  transfer  function  resulting  from  a  stabilizing 
output  feedback  controller  is  derived. 

Fintdly,  since  periodic  controllers  are  known  to  relocate  control  loop  zeros  [7],  we  give  a 
solution  to  the  zero  placement  problem  for  open  and  closed  loop  system  and  characterize  control- 
theoretic  properties  of  the  resulting  system.  In  this  thesis,  we  treat  the  above  problem  by 
considering  a  linear  time-invariant  SISO  system  of  the  form 

x{t)  —  Ax{t)  -{-  Biu{t)  -f  B2w{t)  ,  (0.5) 

y(t)  =  Cx(t)  ,  X  e  nV',u,w^y  &  IR  . 

In  the  design  stage  of  the  system  (0.5),  when  no  actuator  positioning  and  disturbance  protection 
measures  are  yet  finalized,  input  vectors  Bi  and  B2  and  output  vector  C  may  be  viewed  as  free 
parameters  to  be  chosen  so  that  appropriate  performance  specifications  are  satisfied.  Among 
these,  it  seems  reasonable  to  require  that  the  transmission  from  the  control  u  to  the  output  y  be 
maximized  and  the  transmission  from  the  disturbance  u;  to  y  be  minimized.  Since  a  choice  of 
B\  and  B2  defines,  for  a  given  A  and  C,  the  system’s  zeros,  the  problem  of  choosing  B\  and  B2 
also  defines  the  zeros  of  the  system  (0.5).  The  purpose  of  this  research  is  to  give  a  solution  to 
the  zero  placement  problem  for  open  and  closed  loop  system  and  characterize  control-theoretic 
jroperties  of  the  resulting  system.  Note  also  that  zdthough  the  problem  of  optimal  pole  placement 
las  received  enormous  attention  during  the  last  30  years,  that  of  optim^  zero  placement  has 
>een  relatively  neglected  in  control  systems  research.  This  latter  problem  could  however  be  of 
substantial  importance  not  only  in  the  context  of  vibrational  control  but  adso  in  the  placement 
of  actuators  and  sensors  since  the  latter  determines  the  input  £uid  output  vectors  of  the  linear 
model  (0.5)  [8]-[10]. 
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5  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 


5.1  PART  1.  STATE  AND  OUTPUT  FEEDBACK  STABILIZ ABILITY  AND 
POLE  PLACEMENT  CAPABILITIES 

5.1.1  INTRODUCTION 

In  this  part  of  the  report,  we  develop  a  control  theory  for  a  class  of  SISO  dynamical  systems 
described  by  the  following  equations: 

i(t)  =  Ax{t)  +  ,  (1.1) 

y{t)  =  Cx{t)  , 

m  =  /(f  +  r),T^0, 

1 

-jT  /(0d<  =  0, 

0  <  c  <  1  , 

where  A  €  is  a  periodic,  zero  average  scalar  function,  and  e  is 

a  small  positive  constant.  A  characteristic  feature  of  this  problem  is  that  the  control,  u,  enters 
the  open  loop  dynamics  as  amplitude  of  a  periodic,  zero  average  function.  Such  situations 
arise  in  a  number  of  applica.  ons  where  two  conflicting  control  goals  have  to  be  accomplished  by 
a  single  actuator.  For  instance,  in  the  helicopter  control  problem,  a  single  actuator  (the  blades’ 
pitch  angle)  is  used  to  ensure  both  the  desired  altitude  and  the  fuselage  oscillations  suppression. 
These  goals  are  conflicting  in  the  sense  that  if  the  pitch  angle  is  chosen  to  ensure  the  desired 
altitude,  the  fuselage  oscillations  are  not  suppressed;  if  the  pitch  angle  is  used  to  suppress  the 
fuselage  oscillations,  the  desired  altitude  is  not  attained.  In  order  to  accomplish  the  two  goals 
simultaneously,  a  frequency  separation  approach  may  be  employed.  Specifically,  a  low  frequency 
control  may  be  used  to  stabilize  the  desired  altitude  and  a  high  frequency,  average  zero  control 
may  be  used  to  suppress  the  vibrations.  When  the  control  loop  is  closed  with  respect  to  the  low 
frequency  control,  the  resulting  equations  have  the  form  of  system  (1-1),  and  the  goal  is  to  choose 
the  control,  u,  as  a  function  of  x  or  y  so  that  the  system  has  the  desired  dynamical  properties. 
This  ideology  has  been  implemented  in  the  Higher  Harmonic  Control  (HHC)  of  helicopters,  where 
periodic  feathering  of  rotor  blades  around  a  fixed  pitch  angle  is  introduced  in  order  to  suppress 
the  fuselage  vibrations. 

Another  example  is  the  periodic  operation  of  chemical  reactors  [19].  Here  again  the  problem 
is  to  choose  the  amplitude  of  input  flow  vibrations  so  that  the  closed  loop  system  behaves  as 
desired. 

In  order  to  simplify  the  analysis  and  obtain  constructive  results,  following  [11]-[17],  we  assume 
that  the  periodic  function  /(<)  is  of  high  frequency  as  compared  with  the  dynamics  of  i  =  Ax. 
Formally,  this  means  that  function  /  has  the  asymptotic  form  /(^),  where  e  >  0  is  sufficiently 
small.  Thus,  more  precisely  the  problem  addressed  in  this  part  of  the  report  is  as  follows: 

Given  system  (1.1),  determine  under  what  conditions  there  exist  K  €  and 

Co  ■<  1  such  that  for  all  c  <  cq,  the  closed  loop  system  composed  of  (1.1)  with  the  time  invariant 
state  space  controller 

u  =  — X  ,  (1.2) 

c 
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or  a  time  invariant  output  controller 


X  =  Ax  +  Buf  +  L{y  ~  y)  ,  (1.3) 

y  =  Cx  . 

A'  . 

u  =  — X  , 

e 

is  asymptotically  stable.  The  state  feedback  gains  of  (1.2),  (1.3)  are  restricted  to  be  time  invariant 
for  reasons  of  practical  implementation.  Problem  (1.1),  (1.2)  is  considered  in  Section  5.1.2  and 
problem  (1.1),  (1.3)  is  discussed  in  Section  5.1.3.  In  addition,  we  characterize  the  pole  placement 
capabilities  ensured  by  closed  loop  vibrational  control  and  present  the  corresponding  results  in 
Section  5.1.4.  To  illlustrate  the  results,  in  Section  5.1.5  we  consider  an  example  inspired  by  a 
helicopter  with  HHC. 

5.1.2  STATE  SPACE  FEEDBACK 

In  this  section,  we  present  the  result  for  stabilization  of  the  system  (1.1)  with  state  space 
feedback  (1.2). 

Theorem  1.1:  There  exists  a  AT  €  and  an  cq  >  0  such  that  for  all  0  <  e  <  cq 

system  (1.1),  (1.2)  is  asymptotically  stable  if  only  if  (A,  B)  is  stabilizable  and  the  sum  of  all  the 
controllable  eigenvalues  of  A  is  negative. 

Proof:  Necessity  is  proved  by  the  following  considerations.  Represent  the  state  space 
model  (1.1)  in  Kalman  canoniccil  form, 


where  the  pair  {Ac,Bi)  is  controllable.  Since  A„c  is  not  affected  by  feedback,  the  stabilizability 
of  (A,  B)  is  necessary. 

The  stability  of  Xc,  with  x„e(0)  =  0,  is  governed  by  the  following  equation: 

Xe  =  AeXc  + Biuf  ,  Xc  €  IR’",  U  €  IR  .  (1.5) 

Introducing  a  state  feedback  u  =  iifxc/c,  we  obtain 

i.  =  [a,  +  bAj  (^)]  I.  .  (1.6) 

Since  (1.6)  is  periodic,  there  exists  a  Lyapunov  transformation  which  reduces  (1.6)  to  an  equation 
with  constant  coefficients, 

i  =  Az  , 

preserving  the  stability  property.  From  the  Jacobi-Liouville  theorem  [13], 

Tr[A  +  BAf{^]^t=  TrA, 
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where  T  is  the  period  of  Thus, 


Tr  Ac  =  Tr  A  , 

where  Tr  Ac  is  equal  to  the  sum  of  all  the  controllable  eigenvalues.  Therefore.  Tr  Ac  <  0  is 
necessary  and  this  completes  the  proof  of  necessity. 

Suffiency  is  proved  as  follows;  Consider  (1.4)  and  assume  that  all  the  eigenvalues  of  .4ne  ha%'c 
negative  real  parts.  Without  loss  of  generality,  assume  that  (1.5)  is  in  the  controller  canonical 
form,  i.e.; 

Xc  =  AcXc  +  Biuf  ,  Xc  G  JR"',  u  G  JR  , 

where 


0 

1 

•  0 

■  0  ■ 

Ac  = 

0 

0 

•  1 

,  B,= 

0 

-flm 

Urn  — 1 

•  -fll 

.  1  . 

and  Qj  are  the  coefficients  of  the  characteristic  polynomial  of  matrix  Ac-  Apply  state  feedback 


u  = 


(1.7) 


where  A:,  ~  1,  i=2,  . , m.  The  closed  loop  system  is 


Xc  =  AcXc  +  ^BiKf  Xc 

(1.8) 

(1.9) 

The  generating  equation,  [17],  for  this  system  has  the  form 

(1.10) 

where  t  =  tit.  The  general  solution  of  (1.10)  is 

Xc  =  $(t)IcO  , 

where  ^{r)  is  a  fundamental  matrix  for  BK\f{T)  and  ico  is  a  constant.  Consequently,  introducing 
the  substitution 


Xc  =  = 


we  obtain  the  following  equation  in  Bogoliuboff’s  standard  form  [21]; 


dr 


\dh 


^  -1 


$-'(r)/lc$(r)C  . 


(1.11) 
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Applying  the  averaging  principle  [21],  we  obtain  the  following  averaged  equation 


(1.12) 


where 


(D  (D  = 


(1.13) 


Vcy  Vc^”  0  _  0  _  •••  1  ’  '  ' 

.  Om  — 1  ^2^m— '  '  ‘  . 

and  bar  denotes  the  averaged  value,  i.e.  ^(t)  =  1/T  Jq  ^{T)dT.  Let  Ao, , . . . ,  \o„  be  the  open  loop 
eigenvalues  of  (1.5)  and  choose  the  closed-loop  eigenvalues  A^, , . . . ,  A^,,  of  (1.12)  as  follows: 


1  ”* 

Ac^  =  —  52  Aoi  +  j  Im  Ao, 


(1.14) 


Then  the  state  feedback  gains  (1.7)  can  be  found  to  be 

t  ~  ^«)  _  9  „ 

~  I  "72  ’  i  --  2, . . . ,  m  , 


(1.15) 


where  Cc, ,  •  •  •  1 0cm  coeflScients  of  the  closed  loop  characteristic  equation  corresponding  to 

Acj ,  .  .  .  )  Ae„ . 

The  control  gains  (1.15)  guarantee  the  asymptotic  stability  of  the  averaged  system  (1.12).  As 
it  follows  [ITj,  if  (1.12)  is  asymptotically  stable,  there  exists  cq  >  0  such  that  for  all  0  <  c  <  Co 
equation  (1.5)  is  also  asymptotically  stable.  This  proves  the  sufficiency.  Q.E.D. 

Corollary  1.1;  Assume  that  the  sum  of  the  controllable  eigenvalues  of  A  in  (1.1)  is 
positive,  then  there  exists  am  cq  such  that  no  dynamic  state  feedback  of  the  form 


V  =  Fv  +  Gx  , 
u  =  ^  [Hv  +  Jx]  , 


(1.16) 


will  stabilize  the  system  (1.1)  whenever  0  <  c  <  co- 

Proof:  The  resulting  closed  loop  equations  with  the  dynamic  state  feedback  controller 
(1.16)  are  : 

i  =  +  . 

V  =  Fv  -f  Gx  . 


In  fast  time  t  —  if e,  the  above  equation  becomes  : 


I  _  eA  +  BJf( 

S  ”  eG 


r)  BHf{T)  X 
tF  V 


(1.17) 
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Let  $(t)  be  a  fundamental  matrix  for  BJ /{t).  Define 


ViXr)=  r  Hr- q)BHf{q)dq. 
Jo 


and  the  substitution 


[  a:(T) 

*- 

<^(t)  \V{t) 
0  I 


ar) 


Using  the  above  substitution,  we  obtain 


L  dr  i 


=  e 


^-^AW -^-^WGW -^-^WF 

G$  GVU  +  F 


(I.IS) 


which  is  an  equation  in  the  standard  form  [21].  Introduce  the  following  Lyapunov  transformation 


2i(0 

22(t) 

0  / 


m 

q>{t) 


which  will  preserve  the  stability  of  (1.18),  yields 


■  (9  (0  «■'  (0  (i)  ' 

Zl 

.  -^  2  . 

G$  (i)  F 

.  ^2 

Applying  the  averaging  principle  [21],  we  obtain 


$-iA$  0 

G$  F 

[15  J 

Thus  the  eignevalues  of  the  averaged  closed  loop  system  are  the  union  of  those  of  and  F. 

It  follows  from  Theorem  1.1  and  the  averaging  principle  [21],  that  if  Tr  A  >  0,  there  exists  an  cq 
such  that  the  closed  loop  system  (1.1),  (1.16)  will  be  unstable  whenever  0  <  c  <  Cq-  Q.E.D. 

5.1.3  OUTPUT  FEEDBACK 

This  section  deals  with  the  problem  of  stabilizing  system  (1.1)  with  an  output  feedback  (1.3). 
Theorem  1.2:  There  exists  a  K  €  and  an  co  >  0  such  that  for  all 

0  <  c  <  Co  the  system  (1.1),  (1.3)  is  asymptotically  stable  if  and  only  if  (A,  B,G)  is  stabilizable 
and  detectable  and  the  sum  of  the  controllable  eigenv2dues  of  A  is  negative.  The  separation 
principle  holds,  i.e,  the  choice  of  K  and  L  can  be  carried  out  independently. 

Proof:  Necessity  is  proved  by  the  following  considerations.  Consider  the  observer, 

X  =  Ax  +  Buf  +  L{y  -  Cx)  , 


and  the  feedback  law 


K, 

u  =  — I  . 

c 

The  dynamical  equations  for  the  closed  loop  system  are: 
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.ij  [LC  {A- LC)  +  B^f[^)  \[i  \ 

Defining  the  observation  error  e  =  x  —  i.  we  obtain  the  following  equivalent  dynamical  equations: 

e  J  [  0  {A  -  LC)  J  [  e  . 

Thus,  the  stability  of  the  closed  loop  system  depends  on  those  of  A-\-  {\ / e)B K f{t  j e.)  and  A  —  LC . 
Using  the  results  of  Theorem  1.1  completes  the  proof  of  necessity. 

Sufficiency  is  proved  as  follows:  Consider  system  (1.1)  with  {A,B,C)  stabilizable  and  de¬ 
tectable.  In  fast  time  t  =  i/e,  the  resulting  closed-loop  equations  with  output  feedback  (1.3) 
are: 


■  dx  • 

ik 

_ 

.  dr  . 

eA  BKfir) 

eLC  e{A-LC)  + BKfir) 


Let  $(t)  be  a  fundamental  matrix  for  BKfir).  Define 

$i(r)  =  $“^(r)^$(r)  , 

and  the  substitution 


(1.20) 


x(.r)]_\l  BK  £  <l>Mnq)dq 


[  {(’■) 


Using  the  above  substitution  and  following  the  proof  of  Theorem  1.1,  we  reduce  (1,20)  to  the 
standard  form  [21]  and  apply  the  averaging  principle  to  obtain  the  following  averaged  equations: 


^]  _  \  A  +  ^LC  -LC  ^j-A  +  LC-^LC 


$1  - 


■  £ 


(1.21) 


To  simplify  (1-21),  introduce  the  following  transformation 

^ri_f/  0  ]rr 

Z2  I  —I  W  ' 

which  yields 


-{■  A  —  LC  +  ^LC 
0  A-LC 


2l 


(1.22) 


Using  a  construction  similar  to  that  of  (1.10)-(1.12),  one  can  compute  the  state  feedback  gain 
required  to  stabilize  and  assign  all  the  eigenvalues  of  A  —  LC  through  the  choice  of  L  so  that 
the  asymptotic  stability  of  (1.21)  is  guaranteed.  As  it  follows  [17],  if  (^1.21)  is  asymptotically 
stable,  there  exists  an  Cq  >  0  such  that  for  all  0  <  e  <  cq  equation  (1.19)  is  also  asymptotically 
stable.  Q.E.D. 
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5.1.4  POLE  PLACEMENT  CAPABILITIES 


Consider  again  system  (1.1)  with  feedback  (1.2)  and  assume  that 


K 


e 


where  fcj  ~  1,  t  =  1, . . .  ,n.  Thus,  the  closed  loop  system  is 

i  =  +  .  (1.23) 

As  it  is  shown  in  [17]  and  (1-12),  this  equation  can  be  reduced  to  the  averaged  equation, 

t  =  {A  +  B)l,  (1-24) 

where 


x(t)  =  $(t)^(t)  , 

$(r)  is  a  fundamental  matrix  for  BK /(r)  and  r  =  t/e.  As  it  follows  from  [17]  and  [21],  (1.23) 
is  asymptotically  stable  for  sufficiently  small  c  if  the  averaged  system  (1.24)  is  asymptotically 
stable.  Matrix  B  which,  along  with  matrix  A,  defines  the  stability  of  (1.24)  can  be  characterized 
as  follows: 

Theorem  1.3:  Assume  A  and  jB  are  in  the  controller  canonical  form.  Then  there  exists 
Co  such  that  for  all  0  <  c  <  Co,  matrix  B  of  (1.24)  has  the  form: 


B^ 


0  0  ■ 

0  0 
-?2  0  . 


(1.25) 


and 


^2  ki 


where  a,,  s  ~  1,2,...,  are  the  Fourier  coefficients  of  /(t),  i.e. 


OO 

/(’■)  =  +  • 
»=i 


Proof:  Follows  directly  from  Theorem  3  of  [13]. 


Denote  the  characteristic  polynomials  of  A  and  {A  +  B)  of  (1,23)  and  (1.24),  respectively,  by 


po(s)  =  s’' +  a,os”~^  +  oaos”"*  +  •  •  •  +  a„o  ,  (1.26) 

p(s)  =  s’*  +  Cl  s’*"' +  025""’^  H - t-a„.  (1-27) 

It  follows  from  Theorem  1.3  that 

oi  =  flio  j  (1-28) 

^  020  ,  (1-29) 
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Figure  1.1:  Sector  region  Z)((r,u>) 


and  flj,  3  <  j  <  n,  can  be  arbitr^ily  aissigned.  Below  we  analyze  to  what  extent  the  constraints 

(1.28)  and  (1.29)  prevent  the  control  designer  from  assigning  the  closed  loop  eigenvalues  of  the 
averaged  equation  (1.24)  to  a  desired  region  of  the  complex  plane. 

More  specifically,  considering  the  closed  region  D{a,u))  of  Figure  1.1,  our  purpose  is  to  identify 
the  conditions  under  which  we  can  find  n  real  or  complex  numbers  Aj, . . . ,  A„,  occurring  as  pairs 
of  complex  poles,  such  that  they  are  the  roots  of  the  polynomial  (1.27)  which  satisfy  (1.28)  and 

(1.29) .  When  this  is  possible,  we  say  that  (real  or  complex)  pole  assignment  in  the  region  D{<t,u)) 

using  closed  loop  vibrational  control  is  possible.  Note  that  if  the  closed  loop  poles  are  confined 
to  this  region,  then  the  system  modes  of  (1.24),  and  hence  the  averaged  modes  of  (1.23),  will 
have  a  time  constant  smzdler  than  or  equal  lo  —1/<t  and  decay  exponentially  at  a  rate  greater 
than  or  equal  to  — a/v/tr^  +  . 

Since,  closed  loop  vibrational  control  can  only  modify  the  coefficients  a,o,  2  <  t  <  n,  of  (1.26), 
we  will  assume  throughout  this  section  that  n  >  2.  Lemmas  A.1-A.5  are  contained  in  Appendix 
Al. 

Theorem  1.4:  Pole  assignment  in  the  region  D{a,u)  using  closed  loop  vibrational  control 
is  possible,  only  if 

—  aio<na.  (1.30) 


Proof:  If  Ai, . . . ,  A„  €  D{(t,u.')  are  the  roots  of  the  polynomial  (1.27).  then 


~  Qi  —  A,  , 


Re(A.)<<7.  (1.32) 

Equation  (1.30)  is  obtained  by  equating  the  real  parts  in  (1.31)  and  using 
(1-32).  Q.E.D. 

Theorem  1.5:  Real  pole  assignment  in  the  region  Z)((7,  u;)  is  possible  using  dosed  loop 
vibrational  control  feedback  if  and  only  if 


—aio  <  na  , 

n  -  1  2 

^20  <  "2^-  Oio 


(1.33) 

(1.34) 


Proof:  The  necessity  of  (1.33)  follows  from  Theorem  1.4.  The  necessity  of  (1.34)  follows 
from  the  fact  that  when  Ai, . . . ,  A„  are  the  roots  of  (1.27),  the  coefficients  oi  and  oj  are 


i=\ 

<12  =  ^  A,  Aj  . 


(1.35) 


(1.36) 


The  maximum  value  that  (1.36)  can  achieve  subject  to  the  constraint  (1.35),  (1.28)  is  given  by 
Lemma  A.2  and  is  exactly  the  right  hand  side  of  (1.34).  Therefore  it  (1.34)  is  violated,  pole 
assignment  to  the  region  with  real  poles  is  not  possible. 

To  prove  the  sufficiency  of  (1.33),  (1.34),  assume  they  both  hold.  Choose 

\  -  '^2 • 

It  is  immediately  checked  that  these  real  numbers  solve  the  problem  of  pole  assignment  to  the 
region  D(o’,uf),  which  completes  the  proof.  Q.E.D. 

Theorem  1.5  gives  a  simple  solution  of  the  problem  of  pole  assignment  to  the  region  D{a,u) 
with  real  poles.  When  complex  poles  are  allowed,  results  similar  to  Theorem  1.5  are  obtained: 

Theorem  1.6:  Complex  pole  assignment  in  the  region  Z?((t,u))  is  possible  using  closed 
loop  vibrational  control  feedback  if  and  only  if 
(i)  for  n  even 


— «2io  <  , 


(1.37) 


n  <  “10 


((n-l)(T*+ui’), 


(if  |or(  >  u)  ,  (1.38) 


«20  <  (f  +  (1.39) 
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(ii)  for  n  odd  with  one  real  pole  placed  dX  p  <  a 


-flio  <  (n  -  l)o-  +  p  , 

(1.40) 

and 

020  <  2(n  -  l)iT2  +u;) 

(1.41) 

or 

020  <  p(aio  +  p)  +  (a;  cr)^  ^  2 

aio  +  p\ 
a  ) 

(if  |(t1  <  u;)  . 

(1.42) 

Proof:  We  first  consider  the  case  when  n  is  even  and  let  (a|  >  uj.  The  necessity  of  (1.37) 
follows  from  Theorem  1.4.  The  necessity  of  (1.38)  follows  from  the  fact  that  when  Ai, . . . ,  A„  are 
the  roots  of  (1.27),  the  coefficients  oi  and  02  are 

Oi  —  —  1 

t=i 

(1.43) 

02  —  ^  A,Aj  . 

i,}ml 

i>} 

(1.44) 

The  maximum  value  that  (1.44)  can  achieve  subject  to  the  constraint  (1.43),  (1.28)  is  given  by 
Lemma  A.2  and  is  exactly  the  right  hand  side  of  (1.38). 

To  prove  the  sufficiency  of  (1.37),  (1.38),  assume  they  both  hold.  Choose 

AJ  =  (<r  +  j(-l)‘u;)  ,  t  =  l,...,n. 

n<T  '  ' 

It  is  immediately  checked  that  these  complex  numbers  solve  the  problem  of  pole  assignment  to 
the  region  which  completes  the  proof. 

When  n  is  even  and  |<r|  <  w,  the  proof  is  similar  and  follows  immediately  from  Theorem  1.4 
and  Lemma  A.3. 

The  proof  for  the  case  when  n  is  odd  is  completed  by  using  Theorem  1.4,  Lemma  A. 4  (if 
kl  >  ‘*')  and  Lemma  A.5  (if  |<r|  <  w).  Q.E.D. 

Remark  1.1;  If  =  0,  then  conditions  (1.38),  (1.39)  and  (1.41),  (1.42)  in  Theorem 

2.6  are  identical.  Note  also  that  the  expressions  in  the  right  hand  side  of  (1.39)  and  (142) 
have  larger  values  than  the  corresponding  expressions  in  the  right  hand  side  of  (1.38)  and  (1.41) 
respectively.  Consequently,  it  is  easier  to  satisfy  the  condition  for  assigning  complex  poles  into 
the  region  Z?((t,u;)  such  that  <  0  with  closed  loop  vibrational  control. 

Remark  1.2:  Some  observations  concerning  the  region  D(<t,u>)  are  in  order.  Let  A  = 
— Cwn  ±  jud  be  a  complex  pole,  where  C  >8  the  damping  ratio,  0  <  C  ^  l>‘*^n  =  |A|  is  the 
natural  (undamped)  frequency,  and  Ud  —  t*>n\/l  —  is  the  damped  natural  frequency.  Then  if 

A  €  D{<r,uj),  it  follows  that  ^  +  w*  and  — C^n  ^  o’-  To  analyze  the  pole  placement 

capabilities  of  the  closed  loop  vibrational  control,  we  note  that  in  practice,  design  specifications 
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Figure  1.2:  Inverted  pendulum  mounted  on  a  platform 


are  often  given  in  terms  of  >  Cmm  and  — C^n  <  ot-  Such  specifications  will  be  satisfied  by 
Z)(o-,u>)  if  <r  <  a  and  or  equivalently 

<r  <  a  ,  (1.45) 

- 1 .  (1.46) 

Hence,  different  values  of  a  and  u>,  subject  to  the  constraints  in  Theorem  1.6,  can  be  chosen  to 
enforce  different  bounds  on  the  damping  ratio,  natural  frequency  and  damped  natural  frequency 
of  the  averaged  closed  loop  system. 

5.1.5  AN  ILLUSTRATIVE  EXAMPLE 

Below  we  present  an  example  inspired  by  a  helicopter  with  HHC.  The  system  is  shown  in 
Figure  1.2.  Here,  an  inverted  pendulum  is  mounted  on  a  platform,  and  the  goal  is  to  maintain 
the  platform  altitude  at  the  desired  level,  2  =  0,  and  the  pendulum  at  the  upright  position,  0  =  0, 
using  a  single  force  actuator,  F,  located  on  the  platform.  The  two  goals  axe  clearly  conflicting 
in  the  sense  mentioned  in  Section  5.1.1.  Specifically,  if  we  use  the  actuator  to  maintain  the 
platform  at  z  =  0,  the  pendulum  will  topple  down;  if  we  use  the  actuator  to  force  the  platform 
to  move  with  an  acceleration  lesser  than  —g,  the  pendulum  will  stay  in  the  upright  position  but 
the  platform  will  not  be  at  2  =  0.  Therefore,  closed  loop  vibrational  control  is  the  method  of 
choice. 

The  Lagrange  equations  for  the  system  at  hand  are  : 

(mi  +  m2)2  —  m2/0sin0  —  mjW^cos^  +  (mi  +  =  F  —  C,z  ^  (L47) 

m2/*0  —  m2/(2  +  </)sin0  =  —  t;>0  ,  (1-48) 

where  mi,2  and  C  denote  the  mass,  altitude  and  damping  coefficient  of  the  platform,  m2, 0,1, 
and  r)  represent  the  mass,  angle,  length  and  damping  coefficient  of  the  pendulum.  Following  the 
ideology  of  closed  loop  vibrational  control,  we  choose  F  with  frequency-separated  components  : 

F  =  {mi-\-m2)g-kz  +  ^fi(^^  .  (1.49) 
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Here  the  first  term  compensates  in  an  open  loop  fashion  for  the  system's  weight,  the  second  term 
is  the  low  frequency  feedback  to  maintain  r  =  0,  and  the  third  term  is  the  high  frequency  control 
to  maintain  ^  =  0. 

To  reduce  (1.47),  (1.48)  to  the  form  (1.1),  we  linearize  all  the  nonlinearities  in  (1.47)  and 
(1.48)  retaining,  however,  the  interaction  term  zQ  in  (1.48).  Assuming  that  fi(tlt)  =  sin(t/t), 
solving  (1.47)  as  t  oo  and  substituting  the  resulting  periodic  function  in  (1.48),  we  obtain  a 
system  of  the  form  (1.1)  with 


X 


0  1 
i  J 


,B  = 


and 


where 


u  = 


1 


K  = 


—a 

/(m,  +  m2) 


0 


(1.50) 

(1.51) 


and  <f>  is  the  phase  shift  introduced  by  the  system  (1.47)  with  input  f\{t/e\.  Since  Tr  /I  <  0  and 
{A,B]  is  controllable,  according  to  Theorem  1.1,  tnere  exists  an  co  such  tnat  for  all  0  <  e  <  cq, 
the  system  (1.50)  is  stabilizable  by  a  state  space  feedback.  Since,  in  addition,  (C,  A)  is  observable, 
the  system  is  stabilizable  by  output  feedback  as  well  (Theorem  1.2). 

Further  according  to  Theorem  1.5,  real  pole  assignment  in  the  region  D{a,u;)  is  possible  if 
and  only  if 


cr  > 


-7 

2m2/2  ’ 


(1.52) 


The  second  condition  is  always  met.  Therefore,  for  real  pole  assignment,  cr  has  to  be  chosen 
so  that  the  condition  (1.52)  is  satisfied. 

According  to  Theorem  1.6,  complex  pole  assignment  into  the  region  D(<7,u;)  is  possible  if  and 
only  if  (1.52)  is  met  and,  in  addition. 


Zl 

I 


is  satisfied.  Again,  since  the  last  condition  is  always  satisfied,  (1.52)  is  necessary  and  sufficient 
for  both  real  and  complex  pole  assignment  in  D(<t,u). 

To  design  a  specific  controller  for  (1.47),  (1.48),  we  average  equation  (1.50)  to  obtain 


1 


0 

g  _ 

I  2?5(m7+m7F 


1  ■ 


e- 


Then,  if  we  choose  the  closed  loop  poles  at 


A2C  = 


7 

2m2B  ’ 
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the  resulting  control  gain  K  is  given  in  (1.51)  with 


Q  =  \/2/(m,  +  ^  .  (1.53) 

In  general,  if  (1.52)  is  met,  the  gain  K  that  assigns  the  closed  loop  poles  at  Ajc  and  \2c,  is 
expressed  by  (1.51)  with: 


a 


\/2/(mi  +  m2)^XicX 


4-  ^ 

+  y 


(1.54) 


If  we  choose  the  observer  gain  L  =  [/i  /2I  to  ensure  that  the  observer  poles  are  at 


then 


Ajo  —  A20  —  "yX 


Ic 


77 

2m2/2  ’ 


7  >  0  , 


^  _  (7  -  1)7 

^  m2l^  ’ 

,  ^  (7  -  2)^7^  g 

^  4mll^  1  ■ 


In  general,  the  observer  gain  L  that  assigns  the  observer  poles  at  Aio,  A20  is  given  by 


(2 


_A  -JL. 

''lo  ^2o  ^  9 

m2P  /  ■ 


To  illustrate  the  behavior  of  the  pendulum-platform  system  with  closed  loop  vibrational 
control,  we  carried  out  numerical  simulations  of  equations  (1.47),  (1.48)  with  the  following  pa¬ 
rameters  :  mi  =  0.1  kg,  m2  =  0.01  kg,  g  =  9.8  m/s^,  /  =  1  m,  j;  =  0.1,  C  =  0.6.  The  control 
law  has  been  chosen  as  in  (1.49)  with  fc  =  5,a  =  2  (see  (1.54))  amd  c  =  0.01.  The  results  are 
illustrated  in  Figures  1.3  and  1.4  for  the  initial  conditions  ;  =  0  m  and  ^  =  0.1  rad.  In  addition, 
we  introduced  impulsive  perturbations  at  the  platform  (at  t  =  2.5  sec)  and  at  the  pendulum  (at 
t  ~  3.5  sec).  As  it  follows  from  these  graphs,  the  closed  loop  vibrational  control  indeed  ensures, 
on  the  average,  simultaneous  satisfaction  of  the  two  conflicting  objectives  using  a  single  actuator. 
Note  that  the  limit  cycle  in  the  dynamics  of  the  platform  is  due  to  the  fact  that  the  controller 
(1.51),  derived  for  the  linear  system  (1.1),  (1.50),  is,  in  fact,  applied  to  the  original  nonlinear 
system  (1.4‘n,  (1.48).  This  is  why  closed  loop  vibrational  control  may  lead  to  the  satisfaction  of 
the  two  conflicting  go^Lls  not  pointwise  in  time  but  only  on  the  average. 
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Figure  1.3:  Platform  Response  with  Closed  Loop  Vibrational  Control 


Figure  1.4:  Pendulum  Response  with  Closed  Loop  Vibrational  Control 


21 


5.2  PART  2.  STABILITY  ROBUSTNESS  IN  CLOSED  LOOP  VIBRATIONAL 
CONTROL 

5.2.1  INTRODUCTION 

Consider  the  system 

X  =  .  (2.1) 

/(O  =  /(t  +  r),r7^o, 


where  x  €  iR",  u  IR,  and  £  is  a  small  positive  parameter,  0  <  e  <C  1. 

A  control  theory  for  system  (2.1),  referred  to  as  closed  loop  vibrational  control,  has  been 
developed  in  Part  1  of  this  report.  In  particular,  the  following  has  been  proved  : 

Theorem  2.1;  There  exists  K  and  £o  >  0  such  that  for  all  0  <  e  <  £o  system  (2.1)  is 
stabilizable  by  a  state  space  feedback 

tt  =  — z  =  -[A:„  k„_i  ...  fci]z  (2.2) 

if  and  only  if  (A,  B)  is  stabilizable  and  the  sum  of  all  the  controllable  eigenvalues  of  A  is  negative. 

This  section  is  devoted  to  the  property  of  stability  robustness  in  closed  loop  vibrational  control. 
Both  synthesis  and  analysis  problems  are  addressed.  In  the  synthesis  problem,  it  is  assumed  that 
(2.1)  is  the  nominal  plant,  whereas  the  true  plant  is  defined  by 

X  =  (A  +  AA)z  +  Bu/0  ,  (2.3) 

where  AA  is  the  perturbation  matrix.  Assuming  that  the  characteristic  polynomial  of  (2.3) 
belongs  to  a  polytope,  Section  5.2.2  below  gives  a  condition  which  guarantees  the  existence  of 
a  controller  (2.2)  such  that  the  closed  loop  system  (2.3),  (2.2)  is  asymptotically  stable  for  all 
members  of  the  polytope.  In  the  analysis  problem.  Section  5.2.3  gives  a  bound  in  the  spectral 
radius  of  AA  so  that  for  all  AA’s  that  meet  this  bound,  a  controller  (2.2),  which  stabilizes  the 
nominal  plant  (2.1),  also  stabilizes  the  perturbed  system  (2.3).  In  Section  5.2.4,  we  present  a  case 
where  closed  loop  vibrational  control  modifies  the  structure  of  the  perturbation  matrix  AA  and, 
therefore,  may  lead  to  robustness  properties  stronger  than  those  of  conventional  (time-invariant) 
control.  In  Section  5.2.5,  we  will  consider  the  robustness  of  closed  loop  vibrational  control  in  the 
presence  of  high  frequency  unmodeled  dynamics. 

5.2.2  SYNTHESIS 

Assume  that  A  and  B  in  (2.1)  are  in  the  controllable  canonical  form 


■  0  1  ...O' 

■  0  ■ 

A  = 

•  •  • 

0  0  •••  1 

,  B  = 

0 

— flnO  ~a(n-I)0  •  •  •  “Oio 

.  1  . 
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where  a,o  are  the  coefficients  of  the  characteristic  polynomial  of  the  nominal  system  (2.1 ).  .\ssume 
that  the  perturbation  matrix  AA  affects  the  last  row  of  the  system  matrix  .4.  and  that  the 
characteristic  polynomial  of  the  perturbed  system  (2.3)  belongs  to  the  class  V  defined  as  follows; 


V  =  |p(s)  =J2Qjpj{s)  :  Qj  >  0,i  =  J  ‘ 


Pj(s)  =  s"  +  a^s"  ^  +  . . .  +  a^,  j  =  1, . . .  ,m  ,  (2.4) 

are  the  vertex  polynomials  and  aj,i  =  l,...,n  denote  the  ith  coefficient  of  the  ;th  vertex 
polynomial. 

Theorem  2.2:  There  exists  K  and  cq  such  that  for  all  0  <  e  <  cq,  any  perturbed  system 
(2.3)  with  open  loop  characteristic  polynomial  in  V  can  be  stabi.Iized  by  a  single  controller  (2.2) 
if  and  only  if  all  coefficients  a{^j  =  1, . . . ,  m,  defined  in  (2.4),  are  positive. 

The  proof  of  Theorem  2.2,  is  based  on  the  following  lemmas: 

Lemma  2.1:  [61]  Consider  a  polynomial 

p(s)  =  s"  +  +  . . .  +  ,  (2.5) 

where  n  >  3,  >  0,  i  =  1, . . . ,  n.  Let  6*  be  the  positive  real  solution  of  the  equation  5(^+1)^  =  1. 

Then  polynomial  (2.5)  is  Hurwitz  if  the  coefficients  rpi  satisfy  the  condition 


/  =  —  2  . 


Inequality  (2.6)  will  henceforth  be  referred  to  as  the  Lipatov’s  conditon  [61]. 

Lemma  2.2:  Consider  the  following  four  polynomials; 

Pki  (s)  =  5"  +  Ti" +  yt +  73  +  74  5"“'*  +  •  •  • 

Pkiis)  =  s”  +  7r 5""*  +  7f  +  73  +  yt 5”"“  +  •  •  • 

Pk3is)  =  s"+7r5""' +7^5’‘"^  +  73  5""^  +  74"5"~‘‘  +  ... 

Pk4{s)  =  s"  +  7+ s""*  +  7f  s""*  +  73  s""^  +  7^ s^~*  +  . . .  , 

where  n  >  3, 7,“  <  7,^,  1  <  *  <  n  and  7f  >  0.  Let  coefficients  of,  1  <  i  <  n  be  defined  as 


{_  ,  (n-iKi-A)  ,  Li-if 
yx  k  4+4- 

,_in=iX!±21+!li=2i 

k  4+4-7. 


,  if  t  is  odd  , 

,  if  I  is  even  , . 


where  fc  >  0.  Then  there  exists  a  k*  such  that  for  all  0  <  fc  <  fc*,  the  coefficient  aj  is  positive 
and  the  polynomials 


Pkcix  (5)  =  5"  +  (7^  +  Sr(A:))s’‘"^  +  (7^  +  ai{k))s^-'^  +  (73  +  ai(*:))s""^  +  . . .  (2.9) 

Pkci2{s)  =  s”  +  (71"  +^(fc))s'‘"^  +  (7^  +  ai(fc))s"“*  +  (7^  +  a5(fc))s’*“^  +  . . . 

Pket3{s)  =  s"  +  (7,“  +  aT(A:))s"~*  +  (7^  +  aj(A:))5"“^  +  (7^  +  a5(fc))s""^  +  . . . 

PkMs)  =  s"  +  (7i+  +  at(fc))s""*  +  (7f  +  oi(fc))s’‘"^  +  (73  +  ai(fc))s"-^  +  . . . 

are  Hurwitz. 
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Proof:  Let  7^  denotes  the  ith  coefficient,  1  <  z  <  n,  of  the  ;th  polynomial.  1  <  i  <  4.  in 
(2.7).  First,  we  show  that  there  exists  k’  such  that  for  all  0  <  ^'  <  k",  polynomials  (2.9)  satisfy 
Lipatov’s  condition,  i.e. 


(7/-!  +  a/_i(A:)){7/+2  +  ai+2(^')) 

mstx  .  “ 

i<Kn-2.i<><4  (yf  +  at(k))ijf^i  +  ai^iik)) 


<  6’  , 


(2.10) 


where  6*  is  defined  in  Lemma  3.1.  Then,  the  statement  of  this  lemma  follows  directly  from 
Lemma  3.1.  To  show  that  Lipatov’s  condition  holds,  consider 


i<:<4  (7/+^(fc))(7/+r+^(fc)) 


1  <  /  <  n  -  2  . 


(2.11) 


It  is  easy  to  see  that  the  function  fi{k)  is  monotonically  increasing  with  respect  to  k.  Indeed, 
define  ^(k)  =  0  and  70  =  7^  =  1.  Let  c^{k)  =  7+  +a:^{k),c~{k)  =  7"  +  ^(fc),  and  Ac,  = 
-  c~(fc)  =  7+  -  7-  for  0  <  ?  <  n.  Then 


//(fc)  = 


1  < /<n-2  , 


crik)cUk)  ’  -  - 

^  cr-l(^)g<+2(^)  (.  ,  Ac;,i  ^  / 

crwcr+t(^)  V  v 


(2.12) 


Ac/+2 

C,- 2(A^) 


= 


0  cUk))  ( 


1  + 


(2.13) 


where 


■{ 


k^  ,  if  /  =  1  , 
k  ,  if2</<n  —  2. 


(2.14) 


We  will  show  that  each  of  the  three  factors  in  the  right  hand  side  of  (2.13)  is  a  non-decreasing 
function  of  k.  Indeed,  c^  =  IjCf  =  7f  >0,  and  since  c~{k)  =  7,“  -f-  ai(k)  for  2  <  i  <  n,  we 
observe  from  (2.8)  that  the  exponent  of  k  in  cf  {k)  is  negative  or  zero.  Also  Ac,  >  0  for  0  <  9  <  n 
and  from  (2.14),  V’(^i)  <  ^’(^2)  for  0  <  ki  <  k2.  Thus,  fi{k)  is  a  continuous  and  monotonically 
increasing  function  01  k.  Also,  as  it  follows  from  (2.13),  fi{0)  =  0,  and  lim*_oo//(^)  =  00. 
Therefore,  by  the  intermediate  value  theorem,  there  exists  k,  such  that 


MK)  =  ,  1  <  /  <  n  -  2 


(2.15) 


Let  fc*  =  minik*.  Since  fi{k)  is  monotonically  increasing  with  k,  it  follows  that  fi{k*)  < 
//(^r)»  1  ^  ^  ^  ~  2.  Therefore, 

+  r/,.,  ^  C. 

max//(A:)  =  max  — -j  - <  max/, (A  )  <  <5  , 


is  satisfied  for  all  0  <  fc  <  fc*.  Finally,  to  ensure  that  oj  >  0,  i.e. 


fc*  —  72  >  0  , 


(2.16) 


choose  k  in  (2.8)  as  follows: 
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Q.E.D. 


(i)  if  7^  <  0  ,  let  0  <  A:  <  k’  , 

(ii)  if  >  0  ,  let  0  <  ^'  <  min[(7f  ^'*]  . 

Proof:  (of  Theorem  2.2)  Consider  the  closed  loop  system  (2.3),  (2.2)  : 

X  =  {A  +  AA  +  BKf  (J^)x  .  (2.17) 

In  the  fast  time  r  =  t/c,  the  closed  loop  system  (2.17)  is 

^  =  {t{A  +  AA)  +  BKf{T))x.  (2.18) 

at 

Let  $(t)  be  a  fundamental  matrix  for  BK /(t).  Reducing  (2.18)  into  the  standard  form  [21]  and 
then  applying  the  averaging  principle,  we  have  the  following  averaged  system 

j  =  $-1  0  (A  +  A/l)$  0  X,  (2.19) 

where  the  bar  denotes  time  averaging  operation.  According  to  Theorem  2  of  [17],  there  exists 
£o  such  that  for  all  0  <  c  <  co,  system  (2.17)  is  asymptotically  stable  if  the  averaged  system 
(2.19)  is  also  asymptotically  stable.  Moreover,  it  was  shown  in  [52]  that  the  averaged  closed  loop 
system  (2.19)  has  the  following  characteristic  polynomial: 

Pc/(s)  =  s"  +  (oi  +  01)3”“^  +  (oj  4-  02)3"”^  +  .  .  .  +  (Un  +  ^)  1  (2.20) 

where 

37  =  k2ki<l>^  ,  *  =  l,...,n,  (2.21) 

=  0, 

4>{i)  =  [  /{r)  dr  , 

and  the  coefficients  a,-  are  the  coefficients  of  the  characteristic  polynomial  of  the  perturbed  system 

The  necessity  part  is  obvious;  Since  the  ai  coefficient  in  (2.20)  cannot  be  adjusted,  it  is 

necessary  that  a\  >  0,j  =  1, . . .  ,m  for  (2.20)  to  be  Hurwitz. 

Sufficiency  is  based  on  Kharitonov’s  theorem  [58]  and  Lemma  3.2.  For  0  <  n  <  2,  it  is  easy 
to  construct  57, 1  <  i  <  n  to  stabilize  (2.20).  lienee,  we  will  consider  the  case  when  n  >  3. 

Let  7,”  =  minjai,7,^  =  maxjoj,  where  a\  denotes  the  ith  coefficient,  1  <  i  <  n,  of  the  jth 
vertex  polynomial  1  <  J  <  m.  We  first  construct  four  interval  polynomials  which  contain  the 
polynomial  (2.20): 

Pkcllis)  =  3’‘  +  cJ‘3’*-*  +cJs’‘"^+Cj3"”®  +  C4S’‘"‘‘  +  ... 

Pkel2is)  =  S”  +  Cf  3””^  +  Cj 3"“^  +  <4 5”"^  +  c| 3""“  +  .  .  . 

Pkclsis)  =  S"  +  q  3’’“*  +  cj 3”~^  +  4 5""^  +  3"““  +  •  •  • 

PkeiA  (5)  =  3"  +  cj’3""^  +  C7 3”“^  +  C3  3""^  +  4 5”"“*  +  •  •  •  . 

(2.22) 

where 

c~  =  7~  -b  57 ,  1  <  t  <  n  ,  (2.23) 

4  =  4  • 
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Next,  we  need  to  determine  a7  in  (2.23)  satisfying  of  =  O.o^  >  0  so  that  the  four  interval 
polynomials  (2.22)  are  Hurwitz.  It  follows  from  Kharitonov's  theorem  that  if  the  polynomials 
(2.22)  are  Hurwitz,  then  (2.20)  is  also  Hurwitz.  Choose  1  <  z  <  n  as  defined  in  (^2.S). 
then  from  Lemma  3.2,  there  exists  a  k’  such  that  for  all  0  <  k  <  k‘ ,  polynomials  (2.22)  are 
Hurwitz.  The  corresponding  stabilizing  state  feedback  gain  K  can  be  computed  from  (2.21). 
As  it  follows  [17],  for  each  Hurwitz  averaged  closed  loop  polynomial  Pci{s)  in  (2.20)  there  exists 
e-p  >  0  such  that  for  all  0  <  e  <  cr>  the  corresponding  closed  loop  system  of  (2.3),  (2.2)  is  also 
asymptotically  stable.  In  [17],  a  lower  bound  of  cp  was  derived.  This  bound  for  ep  is  a  continuous 
function  of  the  coefficients  of  the  open  loop  characteristic  polynomial  p(s).  Since  the  set  of  open 
loop  characteristic  polynomial  p(s)  E  V  \s  closed  and  bounded,  it  follows  from  the  property  of 
continous  functions  that  a  uniform  lower  bound  of  ep  exists.  The  proof  is  completed  by  setting 
Co  =  minT>  ep.  Q.E.D. 

Remark  2.1  :  The  condition  that  the  coefficient  a{  >  0,j  =  l,...,m,  in  Theorem  2.2  is 
equivalent  to  the  requirement  that  the  trace  of  the  perturbed  matrix  A  +  AA  be  negative. 

Remark  2.2  :  Theorem  2.2  is  an  extension  of  the  result  obtained  in  [62]  for  interval  poly¬ 
nomials.  The  assumption  that  V  is  polytopic  is  weaker  than  Kharitonov’s  interval  polynomial 
assumption  because  it  allows  for  linearly  dependent  coefficient  perturbations. 

Example  2.1  :  Consider  a  6th  order  system  (2.3)  with  open  loop  characteristic  polynomial 
p(s)  €  V  where  'P  is  a  poly  tope  of  polynomials  (convex  hull)  with  the  following  four  vertex 
polynomials 


Pi{s)  =  s®  +  0.5s®  —  0.6s'*  +  1.5s^  +  2.5s^ -h  3.4s  —  3  , 

P2(s)  =  s®  +  0.7s®  -  s*  +  2s^  +  2s^  +  4s  -  1  , 

P3(s)  =  S®  -j-  S®  +  3s^  -f-  3s  -)- 1  , 

P4(s)  —  3^  "h  0.6s®  -f*  0.2s**  l.ls^  +  2.9s^  +  3.8s  -f-  2  . 

This  uncertain  open  loop  system  is  unstable  since  the  vertex  polynomials  Pi(s)  and  P2(s)  are 
obviously  unstable.  With  vibrational  state  feedback  control  (2.2),  the  characteristic  polynomial 
(2.20)  of  the  resulting  averaged  closed  loop  system  is  bounded  by  the  interval  polynomials  (2.22). 
Next,  we  construct  57, 1  <  *  <  6,  as  defined  in  (2.8)  and  determine  k  >  0  so  that  the  interval 
polynomials  (2.22)  are  Hurwitz.  Following  (2.12),  we  define  the  function  fi{k)  as; 


and  hence 


f2{k) 

fsik) 

U{k) 


=  jfc  -f  ^  j  (l  -f  Ac4it®  ®)  , 

=  A:  (l  -f  AC2A:®)  j  , 

=  k(^l  +  j  (1  +  Acek*)  , 


where  Ac,  =  cf  —  c,”  =  max^a^  —  minjC^.  Consequently,  with  7f  =  min^  Oj  =  0.5,6*  =  0.4656, 
we  compute  k*,  1  <  /  <  4,  such  that 


Mk;)  =  {k;rHi+my-')=^\ 


(2.24) 
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(2.25) 

(2.26) 
(2.27) 


Mk;)  =  i-;(i  +  i)(i  +  (fc;)^-*)  =  <5* , 

Mki)  =  k;{\  +  2{k;)^)ii  +  2ki)  =  6- , 

MK)  =  k:{i  +  2{k:y-^){\  +  5{kir)  =  6- . 

Solving  (2.24)-(2.27),  we  obtain  k;  =  0.6536,  JtJ  =  0.2327,^3*  =  0.2925,  k;  =  0.3232.  Thus  k’  = 
min/  k;  =  k^  =  0.2327.  Since  minj  <  0,  for  0  <  A:  <  0.2327,  the  interval  polynomials  (2.22) 
and  consequently,  the  characteristic  polynomial  (2.20)  will  be  stabilized.  Arbitrarily  choose 
k  =  0.23,  then  oi  =  1554.7, oi  =  3.5,5;  =  3237.6,5;  =  -0.8,  and  5;  =  360.3.  With  /(r)  =  sin  r, 
the  corresponding  feedback  gains  in  (2.2)  was  ki  =  0,  fcj  =  55.8,  ^3  =  0.1255,  ^4  =  116.1,^5  = 
—0.0287,  and  k^  =  12.9.  The  stability  of  the  characteristic  polynomial  (2.20)  of  the  averaged 
closed  loop  system  with  the  feedback  K,  can  be  verified  by  performing  eigenvalue  tests  on 
appropriate  Hurwitz  matrices  corresponding  to  each  of  the  vertices  of  the  polytope  of  polynomials 
V  [63].  From  Theorem  2  of  fl7],  the  asymptotic  stability  of  the  averaged  closed  loop  system  will 
ensure  the  aisymptotic  stability  of  the  original  system  (2.3). 


5.2.3  ANALYSIS 

Consider  system  (2.3)  with  feedback  (2.2): 

X  =  {A  +  AA  +  B^f(^^)  X  .  (2.28) 

Assume  that  K  and  e  axe  chosen  according  to  Theorem  2.1  so  that  (2.28)  is  asymptotically 
stable  when  AA  =  0.  Under  this  condition,  we  obtain  the  following  results: 

Theorem  2.3:  There  exists  co  such  that  for  all  0  <  c  <  eo»  the  system  (2.28)  is  aymptot- 
ically  stable  if 

(i)AA*(i))  <  (2.29) 

where  P  is  the  unique  positive  definite  solution  of  the  Lyapunov  equation 

P'AA~^PA2Q  =  0,  (2.30) 

Q  is  some  positive  definite  matrix  and 

A  =  0  0  .  (2.31) 

Proof:  From  (2.19),  the  averaged  closed  loop  system’s  equation  of  (2.28)  is: 


5  =  $-i0(A  +  AA)$0  X, 

=  (A  +  A:4)x,  (2.32) 

where 

SA  =  ♦-.(i)A^*(l). 


and  the  bar  denotes  time  averaging  operation. 
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Note  that  (2.32)  is  obtained  by  the  linearity  of  the  averaging  operation.  From  [07].  we  know 
that  the  averaged  linear  time-invariant  system  (2.32)  is  stable  if 

^max(X4)  <  •  (2-33) 

where  P  is  the  unique  matrix  that  satisfies  the  Lyapunov  equation  (2.30),  Q  is  any  positive 
definite  matrix  and  <Jjnax{^)  represents  the  largest  singular  value  of  M. 

Hence,  for  each  Ai4  satisfying  (2.29),  the  averaged  system  (2.32)  is  asymptotically  stable, 
and  there  exists  such  that  for  all  0  <  c  <  system  (2.28)  is  also  asymptotically  stable. 
In  [17],  it  was  shown  that  is  a  continuous  function  of  the  elements  of  the  matrix  AA.  Since 
the  set  of  AA  satisfying  (2.29)  is  closed  and  bounded,  it  follows  that  there  exists  a  uniform  lower 
bound  for  The  proof  is  completed  by  setting  co  =  min^yi  Q.E.D. 

Note  that  the  perturbation  bounds  given  by  (2.29)  are  indirectly  imposed  on  A.4  through 
the  averaged  perturbed  matrix  AA.  Assume  that  A  and  B  are  in  the  controller  canonical  form, 
then 


where 


AA  =  AA  +  AAi  , 


0 


0  0  I 


AAi 


0  _ 

AAinK-i+i  <f>^ 


0  _  0 


(2.34) 


(2.35) 


and  AAij  denotes  the  (i,  j)th  element  of  the  perturbation  matrix  AA. 

Corollary  2.1:  Assume  that  A  and  B  are  in  controller  canonical  form.  Then  there  exists 
an  eo  such  that  for  all  0  <  c  <  cq,  the  system  (2.28)  is  stable  if 


^max  (A  A)  +  <7n»aa:(^Ai)  < 


^min(Q) 

^max(B) 


(2.36) 


Proof:  The  proof  follows  directly  from  Theorem  2.1,  (2.34)  and  the  fact  that  (Tmax(AA  -|- 

AA^  <  «rTOai(AA)  -F  <T,nai(AAi).  Q.E.D. 

Example  2.2  :  Consider  the  system  (2.28)  with 


A  = 


fo 

1  ■ 

.5  = 

■  0  ■ 

o 

-3 

1 

,/(t)  =  sin(T)  , 


and  K  =  [2  0].  Applying  the  averaging  principle,  we  obtain  the  averaged  system  (2.32)  with 


Choose  Q  =  I,  From  (2.36),  the  allowable  range  of  the  perturbation  matrix  which  ensures  the 
asymptotic  stability  of  the  original  system  (2.28)  is  given  by 

a,„„(AA)-f  2|AAi2|  <  0.382. 
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5.2.4  A  SPECIAL  CASE 

In  a  special  case,  it  is  possible  to  show  that  the  structure  of  AA  may  be  differcnl  from  that 
of  AA.  This  difference  can  be  exploited  in  order  to  derive  stronger  robustness  results  for  closed 
loop  vibrational  control  than  linear  time-invariant  control.  Below  we  give  a  necessary  condition 
that  describes  this  situation. 

Theorem  2.4:  Assume  that  A  and  B  in  (2.1)  are  in  controller  canonical  form.  Then 
vibrational  feedback  (2.2)  modifies  the  structure  of  the  perturbation  matrix  A  A  affecting  the 
averaged  system  (2.32)  only  if  the  first  n  —  1  elements  in  the  last  column  of  AA  are  not  all  zero. 

Proof:  From  equation  (2.35),  we  see  that  the  first  n  —  1  elements  in  the  last  row  of  the 
averaged  perturbation  matrix  AA  can  be  modified  by  the  feedback  gains  i  =  2, . . .  ,n,  only 
when  AAjn  0,  j  =  1, . . . ,  n  —  1.  Q.E.D. 

The  above  result  shows  that  vibrational  control  is  capable  of  changing  some  of  the  elements 
of  the  perturbation  matrix  AA  which  affects  the  averaged  system.  Since  the  asymptotic  stability 
of  the  ave’^aged  system  (2.32)  ensures  the  asymptotic  stability  of  the  original  system  (2.28),  the 
modification  in  the  perturbation  matrix  affecting  the  averaged  system  may  change  the  stability 
robustness  of  the  original  system. 

Indeed,  given  a  Hurwitz  matrix  Aq  and  a  matrix  perturbation  direction  Ai,  it  is  shown  in 
[69]  that  the  largest  interval  (rmim^max)  such  that 


is  strictly  stable  for  all  r  €  (r^ 


■^r  =  Ao  +  r  Ai 
c)  is  given  by 


A-,„(-r(Ao)-‘r(A,))  ’ 


A;i»„x(-r(Ao)-»r(Aj))  ’ 


(2.37) 


where  7’(-)  is  any  linear  mapping  which  transforms  the  stability  problem  into  a  nonsingularity 
problem,  A+„^(M)  denotes  the  maximum  positive  (real)  eigenvalue  of  a  square  matrix  M,  and 
A“,_  (A/)  denotes  the  minimum  negative  (real)  eigenvalue  of  M. 

Since  the  perturbation  bounds  in  (2.37)  depend  on  the  perturbation  direction  Ai,  the  mod¬ 
ification  of  AA  into  AA  in  the  averaged  system  by  closed  loop  vibrational  control  may  lead  to 
improved  stability  robustness.  This  is  illustrated  in  the  following  example  : 

Example  2.3  :  Consider  the  system  (2.3)  with  linear  time-invariant  state  feedback 


A  +  BK 

0  1 
0  0 
.  -6  -11 

■  0  1 
0  0 
-8  -12 


(2.38) 


0  ]  r  0 1 
1  +  0  [  -2 
-6  1 


-1  0 


where  Aq  represents  the  closed  loop  system  matrix.  Assume  a  unidirectional  perturbation 
AA  of  the  form 


matrix 


■  0  1  0  ■ 

AA  =  rAi  =  r  0  0  1  , 

2  10. 

then  the  maximal  perturbation  bounds  (2.37)  under  which  stability  is  preserved  is  given  by 
r€  (-1,4). 
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Next,  we  consider  the  same  open  loop  system  (2.3)  with  /(r)  =  sin  r  and  introduce  closed 
loop  vibrational  feedback  of  the  form  (2.2).  Choose  K  =  [2\/2  \/2  0]  so  that  the  averaged 
closed  loop  matrix  A  in  (2.32)  is  identical  to  .4o.  This  ensures  that  the  dynamics  of  the  system 
(2.3)  with  closed  loop  vibrational  control  is  similar  to  the  linear  time-invariant  state  feedback 
(2.38).  Then  AA  in  (2.32)  will  have  the  form 

_  ’010 

AA  =  r  001 
0  0  1 

From  (2.37),  we  determine  system  (2.32)  is  stable  if  and  only  if  r  6  (  —  1,8),  which  compares  fa¬ 
vorably  with  the  interval  (  —  1,4)  obtained  for  linear  time-invariant  state  feedback.  Unfortunately, 
this  situation  does  not  always  take  place: 

Example  2.4  :  Consider  the  same  nominal  closed  loop  system  as  in  Example  2.3.  However, 
assume  that  A  A  now  has  the  form  : 

'010' 

AA  =  r  0  0  1  . 

2  0  1, 

Then  with  linear  time-invariant  state  feedback  (2.38),  the  maximal  range  which  ensures  the 
stability  of  the  closed  loop  system  is  r  e  (—1,4). 

With  closed  loop  vibrational  feedback  (2.2),  the  averaged  perturbation  matrix  AA  becomes 

_  ’0  10" 

AA  =  r  0  01  , 

0  -1  1 


and  closed  loop  stability  of  the  averaged  system  is  ensured  if  r  €  (—1,3.6298). 

5.2.5  UNMODELED  DYNAMICS 

Generally,  a  model  of  the  system  to  be  controlled  may  become  less  accurate  at  high  frequencies 
because  of  unknown  or  unmodeled  parasitic  dynamics.  Moreover,  these  parasitic  dynamics  may 
change  with  time  or  other  physical  parameters,  and  so  cannot  be  confidently  modeled. 

In  state  space  plant  descriptions,  the  addition  of  high  frequency  parasitic  dynamics  is  called 
a  singular  perturbation,  because  the  perturbed  plant  has  more  states  than  the  plant.  Consider 
the  following  singularly  perturbed  form  of  (2.1) 


X 


Aux  +  Ai2Z  +  Biuf  , 


AjjX  -f  A22Z  -f  B2 


(2.39) 


where  x  €  IR"  is  the  modeled  state,  z  €  H*"  is  the  unmodeled  high  frequency  state,  «  6  K  is 
the  control,  f{t)  is  a  periodic,  average  zero  scalar  function,  /x  is  a  small  positive  constant  and 
0  <  c  <  1. 

Following  (2.2),  we  synthesize  a  state  feedback  of  the  form 


(2.40) 
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where  only  the  modeled  state  x  is  available  for  feedback. 

The  resulting  closed  loop  system  (2.39),  (2.40)  is  given  by 

i  —  -^11(0^  4"  ^^12(0^  )  (2.41) 

fiz  =  A2iii)x  +  A22it):  , 

where 

■^11(0  =  +  /^-j  ,  (2.42) 

Ai2{t)  =  Ai2  , 

■^21(0  =  -^21  +  B2  —  f  , 

■^22(0  =  -^22  • 

In  this  section,  we  determine  conditions  on  /In,  /I12,  A2i.,A22,  K,  and  fi  which  ensure  the  stability 
of  the  system  (2.41).  To  answer  this  question,  we  depend  strongly  on  the  following  result: 
Theorem  2.5:  [70]  Consider  the  linear  time-varying  system 

X  =  An{t)x  +  Ai2{t)z  ,  (2.43) 

fiZ  =  A2iit)x  -f-  A22{t)z  . 

Suppose  AiAt)  are  continuously  differentiable,  Aij{t)  are  bounded,  for  i,j  =  1, 2,  Re  A(y422(0)  <  0 
and  the  reduced  system 

X  -  [Aii(t)  -  ii2(t)Aj2(t)A2i(t)]ar  ,  (2.44) 

is  uniformly  asymptotically  stable.  Then  there  exists  a  positive  number  fio  such  that  whenever 
fi  belongs  to  (0, /io)i  the  system  (2.43)  is  uniformly  asymptotically  stable. 

Inspecting  (2.41),  it  is  clear  that  Aij{t)  are  continuously  differentiable,  i4ij(t)  axe  bounded 
for  i,j  =  1,2.  Furthermore,  the  reduced  system  of  (2.41)  corresponding  to  the  form  of  equation 
(2.44)  is 

i  —  A12A22*  A21  +  —  A12A22  -^2)— /  j  ®  >  (2.45) 

Hence  the  stability  of  the  singularly  perturbed  system  with  closed  loop  vibrational  control  is 
answered  by  the  following  theorem  : 

Theorem  2.6:  Consider  the  system  (2.39).  Suppose  that 

(All  “  Ai2A22A2i,Bi  —  A12A22B2)  is  stabilizable  zmd  the  sum  of  all  the  controllable  eigen¬ 
values  of  All  —  A12AJ2M21  is  negative.  Then  there  exists  a  K,  and  an  co  >  0  such  that  for  all 
0  <  «  <  Co,  system  (2.45)  is  uniformly  asymptotically  stable.  In  addition,  if  A22  is  Hurwitz, 
there  exists  a  positive  number  po  such  that  whenever  fi  belongs  to  (0,  po),  the  closed  loop  system 
(2.39),  (2.40)  is  uniformly  asymptotically  stable. 

Proof:  The  stabilizability  condition  for  (2.45)  follows  directly  from  the  results  of  Theorem 
2.1.  From  Theorem  2.5,  the  uniform  stability  of  (2.45)  and  the  asymptotic  stability  of  A22  will  in 
turn  ensure  the  existence  of  a  po  >  0  such  that  the  closed  loop  system  (2.39),  (2.40)  is  uniformly 
asymptotically  stable  for  all  0  <  p  <  po-  Q.E.D. 

A  conservative  estimate  of  the  value  fiQ  can  be  computed  as  in  [71]  or  [72].  In  particular,  the 
estimate  obtained  by  Kokotovic,  Khalil  and  O’Reilly  [72],  based  on  the  solutions  of  two  Lyapunov 
equations,  is  shown  in  most  cases  to  be  less  conservative  than  that  obtained  by  Javid  [71]. 
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5.3  PART  3.  DESIGN  OF  VIBRATIONAL  CONTROLLER  FOR  PERFORMANCE 
AND  DISTURBANCE  REJECTION 

5.3.1  INTRODUCTION 

Consider  the  class  of  dynamical  systems  described  by  the  following  equations; 

x(t)  =  Ax{t)  +  Bu{t)f  ,  (3.1) 

Z{t)  =  CiX(t)  +  DiuU{t)  +  D:uW{t)  , 

y{t)  =  C  x{t)  +  1 

fit)  =  fit  +  T),T^0, 

0  <  c  <  1  , 

where  x  €  IR!'  is  the  state,  y  €  JR  is  the  measured  output,  r  6  JR  is  the  regulated  output, 
u  €  JR  is  the  control  input,  ly  6  JR  is  the  exogenous  input,  f{t)  is  a  known  periodic,  average 
zero  scalar  function,  and  c  is  a  small  positive  constant.  A  characteristic  feature  of  this  system 
is  that  the  control  enters  the  open  loop  system  dynamics  as  an  amplitude  of  a  periodic,  zero 
average  function,  and  this  amplitude  can  be  chosen  to  depend  on  the  system’s  states  or,  more 
generally,  output.  An  example  of  such  a  system  is  the  helicopter  with  Higher  Harmonic  Control 
(HHC),  where  periodic  feathering  of  rotor  blades  around  a  fixed  pitch  angle  is  introduced  in 
order  to  suppress  the  fuselage  vibrations  [3]-[5].  Another  example  is  the  periodic  operation  of 
chemical  reactors  [19]  where  the  input  flow  vibrations  are  introduced  so  that  the  closed  loop 
system  behaves  as  desired. 

The  theory  for  the  control  of  system  (3.1),  referred  to  as  closed  loop  vibrational  control,  has 
been  developed  in  Section  5.1  of  this  report.  In  particular,  necessary  and  sufficient  conditions  for 
the  output  feedback  stabilizability  of  the  system  (3.1)  with  the  observer-based  output  controller, 
Knom,  defined  by 

x  =  Ax  +  Buf  -I-  L{y  -  y)  ,  (3.2) 

K . 

u  =  — X  , 

y  =  Cx  , 

have  been  established.  However,  the  class  of  all  stabilizing  output  feedback  controllers  has  not 
been  characterized,  although  such  a  characterization  could  be  quite  useful,  e.g.  to  satisfy  desired 
performance  specification  [6]. 

The  purpose  of  this  section  is  to  present  such  a  characterization  using  the 
parametrization  approach.  The  results  obtained  are  quite  similar  to  the  Youla 
parametrization  [731.  Specifically,  we  show  that  the  averaged  closed  loop  tramsfer  function  re¬ 
sulting  from  a  stabilizing  output  feedback  controller  is  an  affine  function  of  an  arbitr2U'y  stable 
transfer  function. 

This  section  is  organized  as  follows  :  Section  5.3.2  determines  the  set  of  all  stabilizing  con¬ 
trollers  for  system  (3.1).  The  parametrization  of  the  corresponding  averaged  closed  loop  transfer 
function  is  described  in  Section  5.3.3.  To  illustrate  the  results.  Section  5.3.4  presents  an  example 
where  the  parametrization  of  the  averaged  closed  loop  transfer  function  is  used  to  design  a  con¬ 
troller  that  achieves  certain  step  response  specifications  for  the  original,  non-averaged  system. 
The  disturbance  decoupling  capability  of  closed  loop  vibrational  control  is  discussed  in  Section 
5.3.5. 
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5.3.2  PARAMETRIZATION  OF  STABILIZING  OUTPUT  CONTROLLERS 


The  following  result  for  stabilization  of  the  system  (3.1)  with  output  feedback  (3.2)  has  been 
derived  in  Section  5.1: 

Theorem  3.1:  There  exists  K.  L  and  an  0  <  £o  1  such  that  for  all  0  <  e  <  co  system 
(3.1)  is  stabilizable  by  the  output  feedback  (3.2)  if  and  only  if  (.4,5,C)  is  stabilizable  and 
detectable  and  the  sum  of  the  controllable  eigenvalues  of  A  is  negative.  Moreover,  the  separation 
principle  holds,  i.e.  the  choice  of  K  and  L  can  be  carried  out  independently. 

Although  (3.2)  contains  a  time  varying  function,  fitje),  it  nevertheless  can  be  viewed  as  a 
time  invariant  controller  since,  firstly,  f{t/e)  is  a  part  ol  the  control  signal  and,  secondly,  K  and 
L  are  constant  gains.  Therefore,  it  is  natural  to  parametrize  all  stabilizing  controllers  for  (3.1) 
in  the  class  of  rational  transfer  functions. 

To  accomplish  this,  augment  the  closed  loop  system  (3.1),  (3.2)  bv  an  auxilliary  controller 
Kq  as  shown  in  Figure  3.1,  where  Kq  is  defined  by  : 

xq  =  Aqxq  +  Bqe  ,  (3.3) 


e  =  y  ~ Cx , 

and  y  and  x  are  defined  in  (3.1)  and  (3.2)  respectively.  Note  that  Kq  is  a  compensator  charac¬ 
terized  by  the  high  gain  i  : 

QM  =  1  [c«(»/  -  AqT'Bq  +  Dq]  . 

The  state  space  realization  of  the  augmented  controller  {KnomtKq)  is  : 

x'e  =  AeXe  +  B^uf  -f  LeiV  -  C^Xf)  ,  (3.4) 

“  =  -  C,I.)  ,  (3.5) 

where 


Xe  = 


,Ae  = 


A  0 
0  Aq 


,J5e  = 


C.  =  \c  O]  ,K,  =  \K  Cq]  ,L,=  g  . 

Theorem  3.2:  Assume  (3.1)  is  internally  stabilized  by  the  nominal  controller  (3.2).  Then 
it  is  also  internally  stabilized  by  the  augmented  controller  (3.4),  (3.5)  if  and  only  if  Aq  is  stable. 

Proof:  The  necessity  is  proved  by  the  following  consideration.  The  internal  dynamical 
equations  for  the  closed  loop  system  with  the  augmented  controller  (3.2),  (3.3)  are 


•  i  I  r  A+B^C/(i)  B^/({)]ri' 

i  =  +  ^  _  iC  +  (Bf  -  B^C)/ (f)  B^/(f)  i 

.^'q\  BqC  -BqC  Aq  JL  . 
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Figure  3.1:  Plant  (4.1)  with  the  augmented  controller  (/Comi 


Defining  the  observation  error  Ci  =  x  —  z,  we  obtain  the  following  equivalent  dyn".mical 
equations  ; 

r  *  ]  r  A  +  -b£/(i)  +  b^c/(x)  ]  r  X  ■ 

c’l  =  0  A  -  LC  0  ei  . 

.  J  0  BqC  Aq  L  ^<3  . 


Thus  the  stability  of  the  closed  loop  system  depends  on  the  stability  of  A-+  BK/efit/e),  Aq  and 
A~LC. 

Sufficiency  is  proved  as  follows  :  Since  the  stability  of  the  closed  loop  system  does  not  depend 
on  Dq,  we  let  Dq  =  0.  The  equations  for  the  controller  (3.4),  (3.5)  can  be  rewritten  as 


Xe  =  (Ae  +  BeKef  Q)  “  LeCt)x^  + 


Ke 

U  =  - Xe 

e 


In  fast  time  r  =  t/e,  the  resulting  closed  loop  equations  with  the  controller  (3.4),  (3.5)  are 


dx 

4  = 

dr  . 


eA 

tLeC  c{A 


BKJ{t)  If  X  ■ 

e-LeCe)  +  BeA'e/(r)  J[xe  ' 
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Let  $i(t)  be  a  fundamental  matrix  for  BeKef{T).  Define 

$(r}  =  , 


i(t) 

x.{t) 


I  [/o]($,(t)-7)1  r  e(T) 

0  ^>,(r)  ip{T) 


Using  the  transformation  (3.9),  we  reduce  (3.8)  into  the  standard  form  [21],  apply  the  aver¬ 
aging  principle  and  obtain  the  following  averaged  equations: 


I  1  ^  r  A  +  [7  0]{^T'LeC  -  L,C)  [I  0]($  -  Ae  - 


Ae  ~  ^LeCe  +  LfCe)  ^ 


(3.10) 

where  the  bar  denotes  time  averaging  operation.  To  simplify  (3.10),  introduce  the  following 
substitution 


which  yields 


l/01(-?+A,  +  4r'i.C.-i.C.)l  ry. 
.  \A-LC  0  1  !i 

“  -BqC  Ag  J  i  ’ 


(3.11) 


where  $2(r)  is  a  fundamental  matrix  for  BK /(t).  Hence,  if  Aq  is  stable,  the  stability  of  the 

original  closed  loop  system  is  ensured  since  K  and  L  are  chosen  to  stabilize  A$2  and  A  —  LC 
for  all  0  <  e  <  €o  [17].  Q.E.D. 

Remark  3.1  :  Note  that  the  augmented  controller  equations  (3.4),  (3.5)  are  actually  the 
equations  of  an  observer-based  controller  for  the  system  (3.1),  with  its  dynamics  augmented  in 
such  a  way  that  the  the  signal  e  is  uncontrollable  from  v.  This  structure  is  the  same  as  the  state 
space  version  of  the  Youla  parametrization  of  all  stabilizing  controllers  for  linear  time-invariant 
systems  (74). 

5.3.3  PARAMETRIZATION  OF  THE  AVERAGED  CLOSED  LOOP  TRANSFER 
FUNCTIONS 

As  it  has  been  shown  in  the  proof  of  Theorem  3.2,  the  dynamics  of  (3.1)  with  the  augmented 
controller  (3.2),  (3.3),  and  Dq  =  0  are  characterized  in  the  average  by  the  equation  (3.10).  In 
this  section,  we  aerive  an  exmicit  parametrization  of  the  averaged  closed  loop  transfer  function 
from  input  to  to  the  averaged  output,  z,  defined  below.  We  also  show  that  this  averaged  output 
approximates  asymptotically  the  actual  non-averaged  output  z  as  e  approaches  zero. 

Assume  that  A  and  B  in  (3.1)  are  in  the  controller  canonical  form,  (if  {A,B)  is  only  stabi- 
lizable,  we  assume  that  the  controllable  part  of  the  Kalman  decomposition  is  in  the  controller 
form)  and  let 


K  =  [k„ 


k\  ]  , 
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where  /fc,  ~  l,i  =  1,. . .  ,n.  The  averaged  closed  loop  equations  (3.10)  with  Dq  =  0.  reduces  to 


f  =  A^  +  [-BK  k20^BCQp, 


LC 

■  A  -  BK  -  LC 

k20^BCq 

BqC 

-BqC 

Aq 

where 

K  =  [  k2kn<jp‘  k2kn-\(p^  •••  k\0^  0  ]  , 

4>{~)  =  /  f{r)dr. 

€  Jo 

Define  the  averaged  output  z  as  follows  : 

2  =  +  Dzuti  +  Dzu,W, 


(3.12) 

(3.13) 


(3.14) 


where  u  =  [—K  k24>^Cqp.  Hence,  from  (3.10)  and  (3.14),  the  averaged  closed  loop  transfer 
function  from  w  to  z  resulting  from  the  augmented  controller  {Knomi  Kq)  can  be  represented  by 


where 


Z  —  {Til  +  Ti2QT2i)w  , 


Tiiis)  Tt,{s) 
T,iis)  0 


=  Ct{sI  —  At)~^Bt  +  Dj  , 


At 


Bj 


Cj 


Dt 


A  -B/C 
LC  A- BK  -  LC  \  ' 

B^  B 
LDyy,  B  ' 

C  -C  \  ' 

Dz\u 

Dy^  0 


(3.15) 


(3.16) 


and  K  and  L  are  the  state  feedback  and  observer  gains  of  the  nominal  estimated-state  feedback 
controller,  Knomi  of  the  averaged  closed  loop  system.  The  averaged  closed  loop  transfer  function 
is  depicted  in  Figure  3.2  where  the  system  (3.1)  and  the  nominal  controller  Knom  are  combined 
into  the  block  T.  _ 

The  stable  transfer  function  (^(a)  has  the  following  state  space  realization 

Ig  =  y4glQ  -|-  Bqt  ,  (3-17) 

t;  =  k2<p^Cqxq  . 

Expression  (3.15)  shows  that  the  averaged  closed  loop  transfer  function  with  Bg  =  0  is  affine 
in  (5(a),  where  (5(s)  is  any  asymptotically  stable  strictly  proper  transfer  function  with  the  gain 


36 


Figure  3.2;  Parametrized  form  of  averaged  closed  loop  transfer  function 

of  order  1  (compare  with  (3.3)).  This  flexibility  in  choosing  (5(s)  can  be  used  to  yield  a  closed 
loop  system  which  satisfies  certain  design  criteria. 

Remark  3.2  :  For  Dq  ^  0,  the  algebra  is  considerably  more  involved,  but  a  similar  result 
can  be  obtained  which  parametrizes  the  averaged  transfer  function  in  terms  of  an  asymptotically 
stable  transfer  function.  _ 

Next  we  establish  the  correspondence  between  the  averaged  output  2  and  the  actual  non- 

averaged  output  z.  j-  /  a  a  \ 

Theorem  3.3:  Assume  that  =  0  and  the  transfer  function  corresponding  to  {A,  B,  Cz) 

in  (3.1)  has  relative  degree  greater  than  or  equal  to  2.  Then,  if  the  averaged  system  (3.10)  is 
asymptotically  stable,  for  any  ^  >  0  there  exists  Co(^)  such  that  for  all  0  <  c  <  eo,  system  (3.1) 
is  also  asymptotically  stable  and  the  following  inequeuity 

||2(<)-^11  <^  ,  <6[0,oo),  (3.18) 

holds. 

Proof:  Let  A  and  B  be  in  the  controller  canonical  form,  (if  (A,  B)  is  only  stabilizable,  we 
assume  that  the  controllable  part  of  the  Kalman  decomposition  is  in  the  controller  form).  From 
(3.1)  and  (3.14), 

||2(0-^ll  =  \\Czx{t)  +  DzMt)-c.'^-D^Mi)h 

where  ^  is  the  averaged  steady  state  of  system  (3.10).  Substituting  relation  (3.9)  into  the  above 
equation,  we  have 

||^(()_;(i)||  =  ||C«-{)  +  C[*,(|)-/|i(e‘"‘'"”-l)BC«l?ll  .  (3-19) 

where 

S=[0  0  •••  0  l]^,C'=[c„  c„_i  •••  C2  0  ]  ,  $2  Q)  =  exp(RAXt/c)). 

Thus  (3.19)  simplifies  to 

IWO-IWII  =  l|C({-f)ll,_ 

<  l<^ll|{ -?ll  • 

0<t<n 

Consequently,  the  result  (3.18)  follows  directly  from  Theorem  1.3  of  [75].  Q.E.D. 
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5.3.4  DESIGN  EXAMPLE 

In  this  section,  we  will  use  the  parametrization  of  the  averaged  closed  loop  transfer  function 
to  design  a  controller  of  the  form  (3.4),  (3.5)  so  that  the  original  system  meets  certain  step 
response  specifications. 

Consider  now  an  example  of  system  (3.1)  with 

■  0  1  0  ]  r  0  1  [O' 

A  =  0  0  1  ,B=  0  ,B^  =  ^  0  , 

0  0  -10  J  [  1  J  [  0 

C.  =  [  10  -1  0  ]  =  [  0  ]  =  [  0  ]  , 

C  =  [-10  1  0  ]  =  [  1  1  ,/(^)  =  •sm(^)  , 

where  z  represents  the  step  response  output  of  the  system  and  y  is  the  measured  tracking  error 
to  a  step  input. 

With  c  =  0.002,  a  controller  (3.2)  that  stabilizes  (3.1)  can  be  defined  by 

K  =  [  3.7032  3.2404  0  ]  ,  L  =  [  -12.5  -75  -50  f  . 

The  corresponding  transfer  functions  Tii(s),  Tui^),  and  T2\{s)  can  be  computed  from  (3.16)  with 

A'  =  [  6  5.25  0  ]  ,  L  =  l  -12.5  -75  -50  f  . 

Assume  that  the  design  specifications  for  the  averaged  closed  loop  system  are  defined  in  terms 
of  the  overshoot,  undershoot,  and  the  rise  time: 

2^  =  sup  z{t)  —  I  <  0.25  ,  (3.20) 

t>o 

zZ  =  sup  — z{t)  <  0.7  ,  (3-21) 

<>o 

true  =  inf  {T  1 1(0  >  0.8  for  t  >  r}  <  1  .  (3.22) 

The  design  specifications  (3.20)-(3.22)  are  closed  loop  quasiconvex  [76],  since  the  set  of  closed 
loop  transfer  functions  satisying  the  design  specifications  is  quasiconvex.  Thus  the  controller 
design  problem  can  be  solved  via  quaisiconvex  optimization.  Following  the  approach  of  [76], 
we  use  the  Ritz  approximation  of  the  augmented  controller  Kn{^)  that  consists  of  the  nominal 
controller  Knomy  and  a  stable  transfer  function  Q(s)  defined  as  a  linear  combination  of  the  fixed 
transfer  functions  C?i(s), . . . ,  Qjv(s)  : 

^  =  i;a.Qi(5),  AT  =  5. 

i=0 

Vector  a  =  [ai . . .  On]  €  JR^,  has  to  be  determined  so  that  the  the  controller  Ks{^)  ensures  the 
desired  closed  loop  specification.  With  T  =  1  in  (3.22),  the  overshoot,  undershoot,  and  the  rise 
time  specifications  are 

5 

X] 0.5. (t)  <  1.25  ,  1.0<t<10.0,  (3.23) 

1 

5 

Y,  QiSi{t)  >  -0.7  ,  0  <  t  <  1.0  ,  (3.24) 

1 
5 

5]a.s.(f)  >  0.8  ,  1.0<f<10.0,  (3.25) 

1 
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i?9  =  0. 

The  resulting  closed  loop  step  response  of  the  averaged  system  with  the  nominal  controller  Knom 
and  the  augmented  controller  Kn{x)  are  shown  in  Figure  3.3. 

The  closed  loop  response  of  the  original  non-averaged  system  (3.1)  with  the  nominal  controller 
Knom  and  the  augmented  controller  {KnomiQ)  are  shown  in  Figure  3.4.  It  can  be  seen  that  the 
actual  and  averaged  closed  loop  step  responses  are  almost  identical  which  is  in  agreement  with 
the  result  of  Theorem  3.3.  Also,  both  the  actual  and  averaged  closed  loop  step  responses  with 
the  respective  augmented  controllers  satisfy  the  design  criteria. 

5.3.5  DISTURBANCE  REJECTION 

The  question  of  when  the  disturbance  w  can  be  completely  decoupled  by  feedback  control 
from  the  regulated  output  z  in  the  system  (3.1)  led  to  the  development  of  geometric  control 
theory.  The  so-called  disturbance  decoupling  problems  for  linear  time-invariant  systems  have 
been  investigated  extensively  in  the  last  two  decades.  The  problem  is  to  find  a  compensator  such 
that  the  closed  loop  transfer  function  from  disturbance  w  to  desired  output  z  is  equal  to  zero. 
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Figure  3.3:  Step  Response  of  Averaged  System 


Figure  3.4:  Step  Response  of  Actual  System 
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Using  the  concept  of  M,  5)-invariance,  the  disturbance  decoupling  problem  with  state  feedback 
(DPP)  was  solved  in  [77].  The  problem  of  disturbance  decoupling  with  state  feedback  and  the 
extra  requirement  of  internal  stability  (DPPS)  was  solved  in  [78],  [79].  A  detailed  reference  for 
the  above  mentioned  problems  can  be  found  in  [80]. 

In  this  section,  we  study  the  disturbance  decoupling  problem  with  closed  loop  vibrational 
control.  In  particular,  for  the  system  (3.1)  with  =  0,  i.e. 


x{t)  =  Ax{t)-\- Bu{t)f  ,  (3.26) 

z{t)  =  CzX{t)  , 

we  will  establish  conditions  under  which  a  state  feedback  of  the  form 

u  =  — X  ,  (3.27) 


can  be  found  which  decouples  z  from  w.  Similar  to  the  approach  in  Section  5.3.3,  we  will  first 
derive  the  averaged  equation  for  the  closed  loop  system  (3.26),  (3.27).  Let  $(t)  be  a  fundamental 
matrix  for  BK /{t)  and  introduce  the  substitution 

^(r)  =  $(t)x(t)  . 

Then  the  corresponding  averaged  equation  of  the  closed  loop  system  (3.26),  (3.27)  is 


k 


z 


(3.28) 


Assume  that  (A,  B)  in  (3.26)  is  controllable,  then  without  loss  of  generality,  let  (A,  B)  be  in 
the  controller  canonical  form: 


■  0  1  ..-o' 

■  0  ■ 

A  = 

0  0  •••  1 

,  B  = 

0 

—On  —On-l  •  •  •  —Oi 

.  1  . 

C,  = 

Co  Cl  •••  Cn-l  ]  . 

With  state  feedback 


«  =  _x  =  -  h  fci  , 

c  e  I  J 

the  averaged  closed  loop  system  (3.28)  reduces  to 

t  =  (A  -  B'K)i  +  B^w  ,  (3.29) 

z  =  Ci, 

where  TT  is  defined  in  (3.12). 

It  follows  from  Theorem  3.3,  that  if  the  transfer  function  corresponding  to 
(A,  B,  Cz)  has  relative  degree  >  2,  the  actual  regulated  output  z  will  be  arbitrarily  close  to 
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r.  Hence,  our  goal  is  to  determine  K  for  the  averaged  linear  time-invariant  system  (3.29)  so  that 
'z  is  decoupled  from  w.  This  will  in  turn  ensures  that  for  the  original  non-averaged  system  (3.26). 
the  regulated  output  z  will  also  be  decoupled  from  the  disturbance  w. 

To  solve  the  disturbance  decoupling  problem  for  the  averaged  system  (3.29),  we  first  introduce 
the  following  definition  ; 

Definition  3.1:  A  subspace  V  C  -V  =  IR"  is  said  to  be  (A,  S)-invariant  (or  controlled- 
invariant)  if  there  exists  K  :  X  —*  U  such  that  {A  —  BK)V  C  V. 

Tneorem  3.4:  [80]  A  subspace  V  C  IR”  is  controlled-invariant  if  and  only  if  AV  C 
V  +  ImB. 

Corollary  3.1:  Consider  the  system  (3.29)  with  C  0.  Then  V  C  kerC  does  not  contain 
the  vector  e„  =  [  0  ■  •  •  0  1  ]^. 

Proof:  We  shall  prove  by  contradiction.  Let  e,  denotes  the  i  —  th  standard  basis  vector. 
Assume  that  Cn  €  V.  From  Theorem  3.4,  for  u,  =  e„,  there  exists  a  tn,  €  V  and  u,  €  C  =  IR  such 
that 

ACn  =  Wi  +  Bui  . 

Hence  iw,-  =  Cn-i  €  V.  Applying  Theorem  3.3  again  with  =  Cn-i,  we  see  that  e„_2  €  V. 
Repeating  this  procedure,  we  see  that  ei,e2, ...,e„  €  V.  Hence  dim(K)  =  n  and  we  obtain  a 
contradiction.  Q.E.D. 

Definition  3.2:  The  system  (3.29)  is  said  to  be  disturbance  decoupled  by  closed  loop  vibra¬ 
tional  control  if  there  exists  K  of  the  form  (3.12)  such  that 

<  A  —  BK\lmBv,  >  C  kerC  , 

where 

<  A|ImB  >=  ImB  +  AlmB  +  . . .  +  A"-^ImB. 

Theorem  3.5;  Assume  that  the  transfer  function  corresponding  to  {A,B,Cz)  in  system 
(3.29)  has  relative  degree  2.  Then  the  system  (3.29)  can  be  disturbance  decoupled  by  closed 
loop  vibrational  control  if  and  only  if  there  exists  a  controlled-invariant  subspace  V  satisfying 
ImBtu  CVC  ker  C. 

Proof;  The  necessity  is  well  known  and  established  in  [80].  ^  prove  sufficiency,  we  need 
to  show  that  there  exists  a,  K  E  ^  ~  BK)V  C  V  and  V  C  kerC} 

satisfying  (3.12).  Let  ui, . . . ,  be  a  basis  for  V.  From  Corollary  4.1,  we  know  that  Cn  ^  V.  Now 
Cn-i  y  since  c„_2  ^  0.  Applying  Theorem  3.4,  for  each  v,,  z  =  1, . . . ,  ^,  there  exists  a  tw,  €  V 
and  Ui  €U  =  IR  such  that 

Avi  =  tPi  +  Bu,  , 

where  B  :  U  -+  X.  Let  v^+i,...,Vn-2,e„_i,e„  be  a  set  of  vectors  such  that 
Vj, . . . , . . . ,  Un-2,  c„  is  a  basis  for  iR".  Define  K  by  the  equation 

[vi  .  •  .  ...  Vfi^2  ®n— 1  “  [^1  •  •  •  ^^+1  •  *  •  ^n— 2  1  ^n]  ? 

where  the  scalars  u^+i, . . .  ,u„_2  are  arbitrary  and  u„_i  >  0  and  u„  =  0.  Clearly,  K  satisfying 
(3.12)  exists  and  is  computable  since  the  matrix  [ui  ...u^u^+i . . . u„_2e„_ie„]  is  non-singular. 

Q.E.D. 

It  follows  from  Theorem  3.3,  that  if  the  averaged  output  z  is  decoupled  from  w,  then  the 
effects  of  w  will  also  be  decoupled  from  the  actual  output  z.  The  corresponding  state  feedback 
gain  K  for  the  original  non-averaged  system  (3.26)  can  be  easily  computed  from  (3.12). 


5.4  PART  4.  VERY  HIGH  HARMONIC  CONTROL  IN  HELICOPTERS 
5.4.1  INTRODUCTION 


Helicopter  vibration  is  a  long  standing  problem.  From  the  earliest  days  of  rotorcraft  develop¬ 
ment,  oscillatory  motions  of  the  non-rotating  portion  of  the  airframe  have  been  a  matter  of  serious 
concern.  Vibration  affects  adversely  pilot  performance  and  passenger  comfort  and  increcises  oper¬ 
ation  costs  and  maintenance  requirements,  limiting  the  use  of  helicopters  in  commercial  aviation 
and  military  operations. 

Traditionally,  rotor  induced  vibration  reduction  efforts  have  focused  on  altering  blade  and 
fuselage  dynamic  characteristics  or  supplementing  the  rotorcraft  with  vibration  isolators  and 
attenuators.  These  passive  methods  have  met  with  limited  success  [3].  Considerable  research 
has  been  done  on  “active”  techniques,  which  make  direct  use  of  vehicle  control.  One  of  the  most 
widely  used  method  of  active  control  is  the  Higher  Harmonic  Control  (HHC).  HHC  involves 
driving  the  blade  pitch  angle  at  the  frequency  N/rev  {N  is  the  number  of  blades)  in  order 
to  cancel  the  effect  of  vibratory  airload  transmitted  to  the  fuselage.  Recent  experiments,  [3]- 
[5],  have  shown  that  HHC  may  lead  to  an  order  of  magnitude  reduction  in  fuselage  vibrations. 
For  this  method  to  be  effective,  not  only  the  amplitude  but  also  the  phase  of  the  oscillations 
introduced  should  be  chosen  appropriately,  depending  on  the  flight  conditions.  The  purpose  of 
this  research  is  to  show  that  if  the  frequency  of  the  oscillations  introduced  is  much  greater  than 
Af/rev,  the  phase  dependency  disappears.  To  distinguish  this  approach  from  HHC,  we  call  it 
VHHC,  the  Very  High  Harmonic  Control.  This  part  of  the  report  presents  a  simple  model  and  a 
hypothetical  explanation  of  the  efficacy  of  VHHC.  The  development  presented  here  is  based  on 
the  ideas  of  vibrational  control.  Vibrational  control  is  a  method  for  changing  dynamic  properties 
of  systems  by  introducing  a  control  which  enters  the  open  loop  dynamics  as  an  amplitude  of 
a  periodic,  zero  average  function.  The  theory  of  this  control  technique  has  been  developed  in 
Section  5.1.  As  it  has  been  pointed  out  in  Section  5.1,  vibrational  control  is  a  useful  tool  for 
achieving  two  conflicting  control  goals  when  only  one  actuator  is  available.  This  is  exactly  the 
case  in  helicopters  where  the  rotor  blade  pitch  is  used  to  accomplish  both  the  desired  aircraft 
altitude  and  the  fuselage  oscillations  suppression. 

The  application  of  closed  loop  vibrational  control  requires  the  knowledge  of  a  dynamic  model 
of  a  system  under  consideration.  Unfortunately,  due  to  uncertainties  associated  with  dynamics, 
measurement  data  and  model  structure,  the  mathematic2d  modeling  of  helicopters  for  handling 
qualities,  performance,  and  flight  control  is  a  very  difficult  problem.  The  modeling  of  helicopters 
is  further  complicated  by  numerous  energy  sources  which  cause  fuselage  vibrations.  These  en¬ 
ergy  sources  are  (i)  alternating  aerodynamics  forces  acting  on  the  rotors,  (ii)  engine  vibrations, 
and  (iii)  aerodynamics  forces  acting  on  the  fuselage  and  nonrotating  parts  of  the  machine.  The 
rotor  system,  which  transmits  the  vibratory  drloads  to  the  fuselage  through  the  rotor  shaft,  is 
one  of  the  most  significant  contributors  to  helicopter  vibrations  [2^.  The  rotating  blades  create 
vibratory  airloads  containing  all  harmonics  of  the  rotor  rotational  frequency  which  are  passed 
from  the  blades  to  the  pylon  and  then  to  the  cabin  through  complicated  load  paths.  These  loads 
are  felt  as  vibratory  forces  and  moments  whose  frequencies  are  integer  multiples  of  the  blade 
passage  frequency  (number  of  blades  times  rotational  frequency).  No  simple  models  describing 
these  phenomena  are  available.  Therefore,  for  the  purposes  of  this  research,  we  develop  a  sim¬ 
plified  qualitative  model  of  a  helicopter  and  consider  only  the  fixed  airframe  vibrations  excited 
by  the  rotor.  Aerodynamic  effects  are  faked  as  follows:  Vibratory  air  load  generated  by  the 
rotating  blades  is  modeled  as  a  periodic  function  of  the  rotational  frequency.  The  aerodynamic 
perturbations  induced  by  the  rotor  on  the  fuselage  are  modeled  also  as  a  periodic  function  with 
frequency  equal  to  the  blade  passage  frequency.  Finally,  incremental  aerodynamic  forces  gener¬ 
ated  by  VHHC  are  modeled  by  another  periodic  function  with  high  frequency.  For  simplicity,  all 
these  functions  are  chosen  as  sinusoids,  however  this  choice  does  not  £iffect  the  qualitative  results 
of  this  thesis. 

Using  this  model,  this  thesis  gives  a  phenomenological  explanation  of  the  efficiacy  of  VHHC 
in  helicopter  oscillations  suppression.  The  outline  of  this  section  is  as  follows:  In  Section  5.4.2, 
a  simple  model  of  a  helicopter  is  developed.  Sections  5.4.3  and  5.4.4  are  devoted  to  analysis  of 
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Figure  4.1:  Simple  model  of  a  helicopter 

this  model  with  and  without  VHHC  in  hover  and  forward  flight,  respectively.  In  Section  5.4.5, 
the  effects  of  wind  gusts  are  examined. 

5.4.2  QUALITATIVE  MODEL 

Consider  a  simple,  qualitative  model  of  a  helicopter  shown  in  Figure  4.1.  The  platform  with 
mass  mi  and  the  pendulum  with  mass  m2  sire  intended  to  model  the  rotor  and  the  fuselage, 
respectively.  The  rotor  generates  the  lift  force,  F,  and  induces  aerodynamic  perturbations,  Ts, 
on  the  rotational  degree  of  freedom  of  the  fuselage.  Force  F  and  torque  Tg  are  modelled  as 

F  =  Fo  +  aosinwt,  .  . 

To  =  sin  Nu}t,  ^ 

where  Fo  is  the  steady  lift  due  to  the  conventional  rotor  blade  pitch  control,  aosincjt  is  the  peri¬ 
odic  part  of  the  lift  due  to  vibratory  airload,  u;  is  the  rotor  rotational  frequency,  N  is  the  number 
of  rotor  blades,  ao  and  ^  are  parameters  that  model  the  strength  of  the  periodic  components. 
When  VHHC  is  introduced,  F  is  defined  as  follows: 

F  =  Fo  +  ao  sin  wf  +  7  sin  I  ,  0  <  c  <  1  .  (4.2) 
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The  small  parameter  £  is  intended  to  model  the  high  frequency  of  VHHC.  Although  in  general 
the  amplitude  of  blade  feathering  is  small  (i.e.  of  the  order  c),  the  aerodynamic  forces  generated 
by  oscillating  airfoils  are  proportional  not  only  to  the  law  of  oscillations  but  also  to  the  first 
and  second  derivative  of  this  law  [83].  Therefore,  the  amplitude  of  the  \'HHC  generated  force 
is  (4.2)  is  assumed  to  be  Note  that  the  effect  of  differentiation,  mentioned  above,  has  been 
observed  experimentally  in  [29]  where  higher  HHC  frequencies  required  smaller  HHC  amplitudes 
to  suppress  the  fuselage  vibrations. 

The  constant  part  of  the  lift,  Fq,  is  chosen  to  compensate  for  the  system  weight  and  to 
produce  a  control  force  to  ensure  the  desired  helicopter  altitude.  Therefore,  Fq  is  defined  as 


(mi  +  m2)g  -  kx 

ro  = - 

cos  /i 


(4.3) 


where  g  is  the  acceleration  of  gravity,  p  is  the  rotor  tilt  angle  and  k  is  the  gain  of  the  altitude 
controller. 

In  addition  to  force  F  and  torque  Te,  two  more  forces  act  on  the  system,  Fx  and  Fy.  These 
are  intended  to  model  wind  gusts  and  are  assumed  to  be  of  an  impulsive  nature.  A  specific  form 
of  Fx  and  Fy  is  described  in  Section  5.4.5. 

Given  the  system  of  Figure  4.1  and  forces  (4.1),  (4.2),  we  write  the  Lagrange  equation  for  the 
platform  mi  (rotor)  and  for  the  pendulum  (fuselage)  as  follows: 


(mi  +  m2)i  +  m^lB  sin  6  +  cos  6  +  (mi  +  m2)^  =  F  cos  g.  —  Fx  —  (ix  ,  (4.4) 

(mi  +  m2)y  +  m2iB  cos  B  —  m2tB^  sin  0  =  F  sin  g  —  Fy  —  (y  ,  (4.5) 

m2i^B  +  m2^(i sinB  +  y cos B)  +  m2ig sin B  =  Te  +  Fxt sin B  +  Fyt cos B  —  rjB  ,  (4.6) 

where  x,  y  and  ^  denote  the  altitude,  horizontal  distance  and  dampling  coefficient  of  the  platform 
and  B,  £  and  t}  the  angle,  length  and  damping  coefficient  of  the  pendulum. 

The  dynamics  of  this  model  are  analyzed  in  Sections  5.4.3-5.4.5  using  both  a  theoretical  study 
and  numerical  simulations.  For  simulation  purposes,  the  parameters  involved  in  (4.4)*(4.6)  have 

been  chosen  as  follows:  mi  =  0.01kg,  m2  =  0.1kg,  £  =  l/\/l0m,  C  =  0.6,  rj  =  0.1,  fc=5,  ^=9.8 
m/sec*,  ao  =  =  0.1,  iV=3,  u;  =  1  and  £  =  0.01. 


5.4.3  ANALYSIS:  HOVER 

To  model  a  helicopter  in  hover,  we  set  g  =  0.  Assuming  that  Fx  =  Fy  =  0  and  considering 
only  the  vertical  and  angular  motions  of  the  pendulum-platform  system,  we  reduce  equations 
(4.4)-(4.6)  to  the  following  form: 

(mi  -f  m2)x  +  m2£B  sin  B  -|-  m2£B^  cos  B  +  (mi  -|-  m2)g  =  F  —  (x  ,  (4.7) 

m2£^B  +  m2^(x  -f  5^)  sin  B  =  Ts  —  rjB  .  (4.8) 


The  effect  of  VHHC  in  suppressing  the  vertical  acceleration  x  and  the  angular  acceleration  B 
of  the  system  (4.7),  (4.8)  is  analyzed  by  considering  three  different  cases.  In  the  first  case,  the 
system  is  constrained  to  move  only  in  the  vertical  direction  (one  degree  of  freedom),  i.e.  ^  =  0. 
In  the  second  case,  the  system  is  analyzed  with  no  constraint  on  its  angular  motion  (two  degrees 
of  freedom)  but  without  VHHC,  i.e.  a  =  0.  FinsJly,  in  the  third  case  we  investigate  the  system 
performance  with  VHHC,  i.e.  o  ^  0. 

5.4.3. 1  One  Degree  of  FVeedom  without  VHHC 

If  we  assume  that  the  pendulum  is  fixed  rigidly  to  the  platform  Md  therefore  cannot  execute 
zmgular  motion,  equations  (4.7)  and  (4.8)  reduce  to  the  following  single  dynamical  equation: 

(mi  -h  1712)1  -f-  (mi  -b  Tn2)g  =  Fo  +  aosinwt  -  ^x  .  (4.9) 


45 


where 


_ 1 _ 

S  2.69  , 


(4.16) 


and 


rjNuj 

m-iig  —  m2PN'^ijj'^ 

Thus,  the  steady  state  angular  acceleration  is  given  by; 

^jj(0  =  —N^(J^Ail3s\n{Nu)t  +  <ti)  .  (4-17) 

Substituting  (4.15)  into  (4.7),  we  obtain  the  steady  state  solution  for  the  altitude  of  the  platform; 

i„(<)  =  Miaosin{uit  +  9i)  —  {N(jjAi0YM2m2icos{2{Nu}t  +  <ti)  +  fli)  ,  (4.18) 


where 


Mi  = 


\J{k  —  (mi  +  m2)^N'^u}^Y  + 

^  0.267  , 


and 


fli  =  —  tan 


■’(rr 


2Nu}( 


(mi  +  m2)47V2u;* 
Thus,  the  steady  state  vertical  acceleration  of  the  system  is; 


) 


(4.19) 


MiOq sin{u:i  +  v’l)  +  4N^J^M2m2{{Ai^Nu)y  cos{2{Nuft  +  ai)  +  fli)  .  (4.20) 

A  comparison  of  (4.20)  with  (4.12)  shows  that  the  two  degrees  of  freedom  system  exhibits 
vertical  oscillations  with  both  the  rotational  and  the  blade  passage  frequencies.  The  latter 
component  contributes  significantly  to  the  fuselage  vibrations  when  the  helicopter  is  in  hover 
[85].  This  leads  to  a  conclusion  that,  in  the  two  degrees  of  freedom  system,  the  aerodynamic 
perturbations,  Tg,  generated  by  the  rotor,  excite  the  fuselage  vibrations  which  then  are  passed 
on  back  to  the  rotor.  This  positive  feedback  is  responsible  for  the  large  vertical  acceleration  of 
the  system  at  hand. 

This  conclusion  is  supported  by  numerical  simulations  of  (4.13),  (4.14).  Results  are  illustrated 
in  Figures  4.3  and  4.4  (a  =  0),  where  equations  (4.13),  (4.14)  are  solved  under  zero  initial 
conditions.  These  fibres,  which  are  in  a  good  agreement  with  estimates  (4.15)-(4.20),  show  that 
aerodynamic  perturbations,  Tg,  excite  oscillations  of  the  pendulum  with  a  frequency  equal  to  the 
blade  passage  frequency  of  Hz.  As  a  result,  the  pendulum  has  an  angular  acceleration  of 
the  same  frequency  and  an  amplitude  of  about  2.7  rad/s^  as  can  be  seen  in  Figure  4.4(b).  The 
dynamic  excitation  of  the  pendulum  is  subsequently  passed  on  to  the  platform  and  generates 
additional  vertical  motion.  Compared  with  Figure  4.2(a),  Figure  4.3(a)  shows  that  the  altitude 
of  the  platform  is  excited  by  both  the  0.5/7r  Hz  vibratory  airload  generated  by  the  rotor  and  the 
S/tt  Hz  vibration  due  to  the  oscillations  of  the  pendulum.  Thus,  as  shown  in  Figure  4.3(b),  the 
resulting  vertical  acceleration  of  the  platform  also  exhibits  the  effects  of  both  the  rotor  vibratory 
airload  and  the  oscillations  of  the  pendulum.  This  latter  component  accounts  for  the  large 
vertical  acceleration  of  about  0.2  m/s^,  i.e.,  an  order  of  magnitude  increase  in  comparison  with 
the  one  degree  of  freedom  system. 
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5. 4.3.3  Two  Degrees  of  Freedom  with  VHHC 

With  VHHC,  the  vertical  and  angular  accelerations  of  the  system  can  be  obtained  by  substi¬ 
tuting  (4.2)  into  equations  (4.7)  and  (4.S): 


X 


e 


mi  +  m2  —  m-i  sin^  6  \  C 

— —  sin  Ntjji  sin  Oqq  sin  ut  +  — sin  (- 
£  e  \e 


A 

—  sin  0  +  m2g  sin^  6  —  m2^  cos  6  —  kx  —  (x 


mi  -f  m2 

mi  -f  m2  —  m2  sin^  0 

Ci 


(sf 


— - -0^  sin  20  +  — - r  sin  0 

+  m2)  £(mi+m2) 


t{mi  +  m2) 


sin  0  —  §  sin  0 - — 

I  mof? 


m2f^  £(mi  +  m2) 


oo  sinuJ^  sin  0 


- r— sin  sin^  H - ^i3sinA^u;<)  . 

+  m2)c  \tJ  m^f}  ) 


(4.21) 


(4.22) 


The  terms  in  the  right  hand  side  of  (4.21)  and  (4.22)  contain  the  aerodynamic  loads  associated 
with  the  rotor  vibratory  airload,  aerodynamic  perturbations,  and  the  VHHC.  The  VHHC  term  is 
of  much  higher  frequency  than  the  other  two  components.  In  practice,  high  frequency  component 
of  the  vertical  acceleration  is  attenuated  by  passive  damping  in  the  cabin.  It  is  therefore  justifiable 
to  consider  only  the  filtered,  or  averaged  with  respect  to  high  frequency,  vertical  and  angular 
accelerations  of  the  system.  This  is  the  standard  methodology  of  vibrational  control  [11]-  [13], 
where  high  frequency  parametric  oscillations  are  introduced  in  order  to  improve  the  averaged  (or 
filtered)  behavior  of  the  system.  To  obtain  the  averaged  equations  we  first  use  the  generating 
equations  technique  of  [20]  that  reduces  (4.21),  (4.22)  to  the  standard  Bogoliubov’s  form  and 
then  use  the  averaging  principle  of  [21],  to  derive  the  averaged  equation.  For  the  system  at  hand 
the  resulting  averaged  equations  are  (see  Appendix  A2  for  details)  : 


^  I '  1  ,  •  2  j 

— — - — ^  —  sm^  +  mjflfsin  <i>- 

TTli  “T  I7I2  “  TTI2  Sin  Y  \  ^ 

\2f2(mi  -1-  m2  —  m2  sin^  <^)*  / 

cos  <f>  —  kz  —  (^z  ~  ^  sin  Nut  sin  ^  -f  Qo  sinwt^  , 


mi  -f  m2  /  m2 

mi  -b  m2  —  m2  sin^  4>  \2(mi  -f  m2) 


/  — o*  sin^  4> 

\2£*(mi  -f  m2  —  m2 


m2sin^  4>Y 


— - r-j-—  +  sin  2^ 

m2  —  m2  sm  <p)  J 


2m-iP[m\  -f  m2  —  m2  sin*  0)  / 

9  •  1  V  ’1 

+<(m.  ;  m,l  ^  -  7  ~ 


(4.23) 


(4.24) 


where  z  and  <i>  denote  the  averaged  altitude  and  angle,  respectively. 

Comparing  (4.21),  (4.22)  and  (4.23),  (4-2^,  we  oWrve  that  the  introduction  of  VHHC  results 
in  additional  terms  containing  amplitude  a.  To  analyze  the  effects  of  these  terms,  we  first  reduce 


(4.23),  (4.24)  to  the  form  (4.7),  (4.8).  Then,  since  ^  is  small,  we  negiect  the  first  term  containing 
in  (4.24)  and  approximate  the  denominator  of  the  remaining  term  containing  by  mi  -)-  m2. 
Hence  we  write  (4.23),  (4.24)  as: 

(mi  -f  1712)2  +  m2^<^sin  <p  +  1712(6^ cos<i>  =  Qosinu’t  —  kz  —  ^z  ,  (4.25) 

m2i^^  +  {  2  +  5  +  TT.,—  I  sin  ci  =  /3  sin  Nujt  -  rjcp  .  (4.26) 

\  2l(mi4-m2)^/ 

Equation  (4.25)  is  identical  to  (4.71  but  unlike  (4.8),  equation  (4.26)  hits  an  additional  term 
containing  a.  Since  the  vertical  acceleration  z  of  the  fuselage  is  generally  of  the  order  of  O.lg  [84]. 
for  small  <i>  we  can  neglect  the  term  m2izsin4>  and  replace  sin^  by  <f)  and  cos<;^  by  1.  Then  the 
steady  state  solution  of  (4.8)  can  be  represented  as 

5j,(t)  =  A2^s\n{Nujt  +  a'2)  , 

where 

A2  =  -  r  •  - -  •  »  (4-27) 

+  2(T^^  -  m2PN^u^)^  + 

^  0.23,  for  a  =  1  , 


o‘2  =  —  tan* 


r)Nuj 


A  comparison  of  (4.27)  with  (4.16)  shows  that  the  VHHC  term  in  (4.26)  leads  to  an  order 
of  magnitude  reduction  of  the  angular  vibrations  of  the  fuselage,  i.e.  to  the  stiffening  of  the 
system.  As  it  follows  from  (4.20)  (with  Ai,  replaced  by  A2  and  ai  by  <72),  this  results  in  an 
order  of  magnitude  reduction  in  the  vertical  acceleration  of  the  system.  Therefore,  the  effect  of 
VHHC  can  be  interpreted  as  the  decoupling  of  the  positive  feedback  from  the  rotor  through  the 
aerodynamic  perturbations  to  the  fuselage  and  then  back  to  the  rotor. 

This  conclusion  is  supported  by  numerical  investigation  of  (4.23),  (4.241  (see  Figures  4.3  and 
4.4  for  a  ^  0):  Figure  4.4(a)  shows  that  VHHC  significantly  reduces  tne  oscillations  in  the 
pendulum  induced  by  the  rotor  vibratory  airload.  The  resulting  vibratory  load  suppression  in 
the  angular  acceleration  is  illustrated  in  Figure  4.4(b).  The  attenuation  of  the  oscillations  in  the 
pendulum  produces  a  reduction  in  the  dynamic  loaa  transferred  back  to  the  platform.  As  can 
be  seen  in  Figure  4.3(a),  the  3/ir  Hz  vibrations  generated  by  the  oscillations  in  the  pendulum  on 
the  vertical  motion  are  neatly  attenuated.  Note  that  this  attenuation  improves  as  the  amplitude 
of  VHHC  is  increased  (see  Figure  4.5).  For  a  =  1,  the  effect  of  these  oscillations  is  virtually 
eliminated.  In  particular,  the  steady  state  vertical  acceleration  is  about  0.02  m/s^.  The  O/tt 
Hz  vibration  due  to  the  oscillations  in  the  pendulum  has  been  eliminated  and  the  remaining 
vibration  in  the  vertical  acceleration  is  at  a  frequency  of  0.5/x  due  to  the  rotor  vibratory  airload 
aosinwt.  Thus,  the  introduction  of  VHHC  forced  the  model  to  behave  almost  as  a  one  degree 
of  freedom  system. 


5.4.4  ANALYSIS:  FORWARD  FLIGHT 


To  study  the  performance  of  a  VHHC-equipped  helicopter  in  forward  flight,  we  consider 
equations  (4.4)-(4.6)  and  assume  that  there  is  no  wind  gusts,  i.e.  Fx  =  Fy  =  0.  Similar  to  Section 
5.4.3,  we  analyze  the  efficacy  of  VHHC  in  suppressing  the  vertical  acceleration  i,  horizontal 
acceleration  y  and  angular  acceleration  6  by  considering  three  different  cases;  (i)  the  system  is 
constrained  to  move  only  in  the  vertical  and  horizontal  direction  (two  degrees  of  freedom),  (ii)  the 
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system  has  no  constraints  on  its  angular  motion  (three  degrees  of  freedom)  but  without  \'HHC. 
i.e.  Q  =  0,  and  (iii)  the  system  with  three  degrees  of  freedom  and  with  \’HHC,  i.e.  q  7^  0.  We 
assume  that  the  system’s  parameters  are  the  same  as  those  defined  in  Section  5.4.2  and  choose 
/i  =  5  deg. 

5.4.4. 1  Two  Degrees  of  Freedom  without  VHHC 

If  the  system  (4.4)-(4.6),  with  Fx  =  Fy  —  0,  cannot  execute  angular  motion,  then  the  system's 
dynamics  are  governed  by  the  following  two  equations: 

(mi  +  m2)x  +  (mi  -I- 1712)9  =  (Fq  +  ao  sin  a;/)  cos//  —  ,  (4. 28) 

(mi  +  m2)y  =  (Fq  +  oq  sin ujt)  sin  fi  -  .  (4.29) 

Here  again  we  observe  the  absence  of  the  coupling  terms  due  to  angular  motion.  The  response  of 
the  vertical  and  horizontal  acceleration  is  affected  only  by  Fo  and  oq  sin  ut  and  no  blade  passage 
frequency  is  involved.  Numerical  solutions  of  (4.28),  (4.29)  are  illustrated  in  Figures  4.6  and  5.7 
for  0.155  m/sec  initial  horizontal  velocity  and  all  other  initial  conditions  equal  to  zero.  From 
these  figures  we  conclude  that  steady  state  vertical  and  horizonteil  accelerations  are  small  (about 
0.02  m/s*  and  0.002  m/s*  respectively)  and  are  of  0.5/7r  Hz  frequency  due  to  the  rotor  vibratory 
airload  aosinut. 

5. 4. 4. 2  Three  Degrees  of  Freedom  without  VHHC 

Consider  equations  (4.4)-(4.6)  with  angular  motion  but  without  VHHC  (i.e.  a  =  0).  The 
vertical,  horizontal  and  angular  accelerations  can  be  expressed  as: 
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(4.30) 

where  F  is  given  in  (4.1).  An  approximate  analysis  of  these  equations  can  be  carried  out 
by  using  the  method  of  Section  5.4.3.2.  Here,  however,  we  investigate  the  response  of  system 
(4.30)  by  performing  simulation  of  these  nonlinear  equations  with  an  initial  horizontal  velocity 
of  0.155  m/sec  and  ^1  other  initial  conditions  equal  to  zero.  Similar  to  Figure  4.4(a),  in  Figure 
4.10(a)  with  a  =  0  we  observe  that  the  aerodynamic  perturbation  T9  generates  oscillations  in  the 
pendulum  with  a  frequency  of  l.b/n  Hz.  These  oscillation  are,  in  turn,  transferred  back  to  the 
platform  and  produce  additional  motion  in  the  vertical  and  horizontal  directions.  Thus,  as  shown 
in  Fibres  4.8(a)  and  4.9(a),  the  resulting  steady  state  altitude  and  horizontal  velocity  display 
the  effects  of  both  vibratory  airload  aosinwt  and  the  oscillations  of  the  pendulum.  Again,  this 
latter  component  accounts  for  most  of  the  vertical  and  horizontal  accelerations. 

5. 4. 4. 3  Three  Degrees  of  fVeedom  with  VHHC 

With  VHHC,  the  system’s  dynamic  equations  are  governed  by  equations  (4.30)  with  F  de¬ 
fined  in  (4.2).  As  mentioned  in  Section  5.4. 3.3,  the  high  frequency  component  due  to  VHHC 
is  filtered  out  by  passive  damping  of  the  cabin,  and  hence  it  is  of  interest  to  consider  the  fil¬ 
tered  vertical,  horizontal  and  angular  accelerations  of  the  system.  Using  the  generating  equation 
method  proposed  in  [20],  we  reduce  (4.30)  to  the  standard  Bogoliubov  form  and  apply  averaging 
analysis  [21].  The  averaged  equations  obtained  are  as  follows  (see  Appendix  A2  for  details)  : 
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(4.31) 

where  z,  u  and  </>  represent  the  averaged  altitude,  horizontal  distance  and  angle  respectively. 

The  effects  of  VHHC  when  a  helicopter  is  in  forward  flight  and  when  a  helicopter  is  in  hover 
are  very  similar:  VHHC  suppresses  the  oscillations  of  the  pendulum  induced  by  aerodynamic 
perturbations  7®;  this,  in  turn,  leads  to  a  reduction  of  the  dynamic  load  passed  on  back  to  the 
vertical  and  horizontal  motions.  This  can  be  verified  by  performing  the  same  analysis  as  in  Section 
5.4. 3. 3.  The  effectiveness  of  VHHC  in  suppressing  the  vertical  and  horizontal  accelerations  when 
the  helicopter  is  in  forward  flights  is  also  illustrated  by  simulating  nonlinear  equation  (4.31)  with 
the  initial  horizontal  velocity  of  0.155  m/sec  and  all  other  initial  conditions  equal  to  zero.  For 
clarity  of  illustrations,  we  present  here  the  simulations  only  for  a  =  0.25  (other  values  of  a  are 
characterized  in  Figures  4.11  (a)  and  (b)).  As  it  can  be  seen  from  Figures  4.8  and  4.9,  the  Z/ir 
Hz  vibrations  in  the  vertical  and  horizontal  motion  are  greatly  attenuated.  By  comparing  with 
Figures  4.6  and  4.7,  we  notice  that  VHHC  stiffens  the  pendulum  and  forces  the  model  to  behave 
as  a  2  degrees  of  freedom  system.  We  observe  in  Figure  4.10(a)  that  with  VHHC  the  steady  state 
angular  position  of  the  pendulum  is  perpendicular  to  the  platform.  Since  the  angular  position  of 
the  pendulum  is  initially  at  zero,  this  accounts  for  the  large  transient  overshoots  in  Figures  4.8 
2knd  4.9. 

5.4.5  EFFECTS  OF  WIND  GUSTS 

Next  we  analyze  the  effects  of  wind  gusts  on  the  performance  of  VHHC-equipped  helicopters 
during  forward  flight.  We  introduce  impulsive  perturbations  F*  (at  t  :=  1  sec  and  t  =  4  sec) 
and  Fy  (at  t  =  2.5  sec  and  t  =  4  sec)  into  the  system;  the  pulse  duration  and  amplitude, 
respectively,  are:  0.5  sec  and  0.1  N.  As  it  follows  from  Figures  4.12  and  4.13,  VHHC  still  cancels 
the  vibrations.  However,  since  the  gusts  Fx  and  Fy  enter  the  system’s  equations  directly  in  (4.4) 
and  (4.5),  VHHC  is  not  very  effective  in  suppressing  their  effects  on  the  vertical  2ind  horizontal 
accelerations  of  the  system. 
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Acceleration 


Time  (sec) 

Figure  4.3:  (b)  Vertical  acceleration  of  helicopter  in  hover  (2  degrees  of  freedom) 
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Acceleration  (rad/s''2)  Angle  (rad) 


Figure  4.4:  (a)  Angle  of  the  fuselage  in  hover  (2  degrees  of  freedom) 


Figure  4.4:  (b)  Anglular  acceleration  of  the  fuselage  in  hover  (2  degrees  of  freedom) 


54 


Alpha 

Figure  4.5:  Vertical  acceleration  of  helicopter  in  hover  for  different  a 
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Acceleration  (tn/s^2)  Velocity  (m/s) 


Figure  4.7:  (a)  Horizontal  velocity  of  helicopter  in  flight  (2  degrees  of  freedom) 


Figure  4.7;  (b)  Horizontal  acceleration  of  helicopter  in  flight  (2  degrees  of  freedom) 
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Acceleration  (in/s^2) 
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Acceleration  (rad/i^2)  Angle  (rad) 


Acceleration  Acceleration  (tnA''2) 


Figure  4.11:  (a)  Vertical  acceleration  of  helicopter  in  flight  for  different  a 


Figure  4.11:  (b)  Horizontal  acceleration  of  helicopter  in  flight  for  different  q 
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5.5  PART  5.  //2-OPTIMAL  ZEROS 
5.5.1  INTRODUCTION 


Another  important  aispect  in  the  control  design  of  systems,  besides  ensuring  sufficient  number 
of  actuators  as  discussed  in  Section  5.1,  is  to  determine  the  location  of  the  control  system 
components,  namely  actuators  and  sensors.  In  this  part  of  the  report,  we  treat  this  problem  by 
considering  a  linear  time-invariant  SISO  system  of  the  form 

x{t)  =  Ax{i)  +  Biu{t)  +  B2w{t)  ,  (5.1) 

y{t)  =  Cx{t)  ,  X  £  IRJ'^u^w^y  ^  IR. 

In  the  design  stage  of  the  system  (5.1),  when  no  actuator  positioning  and  disturbance  protection 
measures  are  yet  finalized,  input  vectors  Bi  and  B^  may  be  viewed  as  free  parameters  to  be  chosen 
so  that  appropriate  performance  specifications  are  satisfied.  Among  these,  it  seems  reasonable 
to  require  that  the  transmission  from  the  control  u  to  the  output  y  be  maximized  and  the 
trtinsmission  from  the  disturbance  10  to  y  be  minimized.  The  term  “transmission”  could  be 
formalized  in  a  number  of  ways.  We  choose  here  a  formalization  in  terms  of  transfer  functions. 
Specifically,  we  say  that  the  transmission  from  u  to  y  is  maximized  if  the  ff2-norm  of  G„(s)  := 
C{sl  —  A)~^Bi  is  maximized  with  respect  to  all  Bi  €  such  that  |lRi||  =  1  (j|  •  [1  denotes 
the  Eucliaean  norm  of  a  vector).  Similarly,  the  transmission  from  u;  to  y  is  minimized  (i.e.  the 
disturbance  rejection  is  optimized)  if  the  /f2-norm  of  G^{s)  :=  C{sl  —  A)~^B2  is  minimized  with 
respect  to  all  B2  €  JR”  such  that  ||B2il  =  1*  Since  a  choice  of  B\  and  B2  defines,  for  a  given 
A  and  C,  the  system’s  zeros,  optimal  B\  and  B2  define  the  zeros,  optimal  with  regard  to  the 
above  stated  criteria.  The  purpose  of  this  research  is  to  give  a  solution  to  the  i/2*optimal  zero 
placement  problem  for  open  and  closed  loop  system  (Sections  5.5.2  and  5.5.4)  and  characterize 
control-theoretic  properties  of  the  resulting  system  (Section  5.5.3).  In  addition,  we  analyze  the 
relationship  between  the  open  and  closed  loop  /fa-optimal  zeros  (Section  5.5.5). 

Since  the  Xoo  norm  of  y  is  often  used  for  control  system  performance  specifications,  the  H2- 
norm  as  a  measure  of  optimality  seems  quite  reasonable  because  it  is  induced  by  the  L2  norm 
on  the  input  and  the  norm  on  the  output  [86].  Other  measures  of  optimality  are,  however, 
possible.  For  instance,  Hughes  and  Skelton  (87)  optimize  the  controllability  and  observability 
“norms”  to  ensure  good  actuators  and  sensors  configuration.  Optimization  of  controllability, 
observability  or  the  information  matrix  has  been  addressed  in  [39],  [90],  and  [91].  An  LQ-type 
approach  is  developed  in  [48].  In  [38],  actuator  location  has  been  addressed  under  a  control  effort 
saturation  assumption.  The  /f2-optimization  of  this  paper  is  related  to  the  above  cited  literature 
because  sensor  and  actuator  locations  typically  determine  the  matrices  C  and  Bi  respectively.  In 
the  present  work,  however,  we  assume  that  the  directions  of  Bi  and  B2  can  be  chosen  arbitrarily 
whereas  in  prar:tical  applications  components  of  vectors  Bi  and  B2  might  be  functionally  related. 
Accounting  for  these  f^unctional  relationships  will  be  a  topic  of  future  work. 

5.5.2  ^2-OPTIMAL  zeros  IN  OPEN  LOOP  ENVIRONMENT 
5. 5. 2.1  Problem  Formulation 

Consider  an  asymptoticadly  stable  SISO  plant  defined  by  a  transfer  function  G(s)  =  C{sl  — 
A)~^B.  Assume  that  vector  B  is  subject  to  the  constraint 

B'^MB  =  1  ,  (5.2) 

where  M  is  a  symmetric,  positive  definite  matrix. 

The  problem  of  ^2-optinial  zero  placement  in  open  loop  environment  is  formulated  as  follows: 

Find  vector  B,  satisfying  (5.2),  such  that  ||G||2  is  maximized  (minimized). 

Remark  5.1:  A  natural  choice  for  M  in  (5.2)  seems  to  be  A/  =  /,  resulting  in  |j5||  =  1.  It 
turns  out,  however,  that  this  normalization  may  lead  to  a  coordinate  dependent  optimal  B  and 
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consequently  coordinate  dependent  zeros  (see  Section  5. 5. 3. 3).  If  (C. -4)  is  observable,  a  better 
choice  is 

M  =  +  A^C'^CA  +  . . .  +  fC'^CA^-^  . 

resulting  in  the  normalization 

{CBf  +  {CABf  +  . . .  +  =  1  .  (5.3) 

Since  CB,CAB, . .  .,CA^~^B  are  the  Markov  parameters  of  G(s)  and  since 

d'git)\ 


CA'B  = 


dv 


,  i  —  0,1,. ..,n  1  , 


«=o+ 


where 

g{t)  =  C-^{G{s)}, 

the  normalization  (5.3)  constrains  B  to  those  that  satisfy 
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(5.4) 


Thus,  the  optimal  B  is  selected  from  all  those  that  result  in  an  impulse  response  which  at  t  =  0"^ 
satisfy  the  normalization  (5.4).  This  in  fswrt  defines  a  unit  sphere  in  the  space  of  impulse  responses 
of  SISO,  LTI  systems  of  dynamic  order  at  most  n.  It  turns  out  that  this  normalization  gives 
a  coordinate  independent  optimal  B  and  as  a  consequence,  coordinate  independent  zeros  (see 
Section  5.5.3.3). 

5.5. 2.2  Main  Result 

Theorem  5.1:  Assume  that  A  in  (5.1)  is  Hurwitz  and  M  in  (5.2)  is  positive  definite 
and  symmetric.  Let  Lo  be  the  observability  Gramian,  defined  by  AJLo  +  LoA  +  C^C  =  0. 

Then  max  IlGullj  (  min  110,^,112  I  is  attained  at  Bi*  {B2*)  collinear  with  the  eigenvector 
Bf  \BfAfB2=l  f 

of  M~^Lo  corresponding  to  its  largest  (smallest)  eigenvalue  AmaxC^j^in)-  Under  this  choice 

Proof:  First,  note  that  M~^Lo  is  a  product  of  two  symmetric,  positive  definite  matrices 
and  so  its  eigenvalues  are  real  and  positive.  Second,  we  observe  that  since  ||G||2  =  B^LoB  (see 
[92])  and  since  (•)*  is  a  monotonic  operation,  the  Lagrangian  for  the  problem  at  hand  is 


L(R,7,  A)  =  -fiB'^LoB)  -  \{B'^MB  -  1)  , 


where  the  multipliers  'y,X  ^  Ul  are  not  simultaneously  zero.  The  first  order  necessary  condition 
is  [93] 

VbL  =  "fLoB  —  XMB  =  0  , 


which  yields 


7M-*LoB  =  AR  . 


Since  B^MB  =  1,7  =  0  implies  A  =  0.  Hence,  it  can  be  assumed  without  loss  of  generality 
that  7  =  1.  Furthermore,  A  >  0.  Thus,  any  optimal  B*  should  lie  along  the  direction  of  an 
eigenvector  of  M~^Lo  and 


||G*||^  =  B-'^LoB-  =  XB-^MB‘  =  X  . 
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Hence, 


max 
B^MB,  =1 

min 


IIG.II^ 


BfLoB;  =  XmaxB-jMBl  =  Amax  , 
BfL.B;  =  X^-,^BfMB;  =  A^i„  . 


Q.E.D. 

Remark  5.2:  Theorem  5.1  can  be  reformulated  for  choosing  the  measurement  vector  C 
to  maximize  (minimize)  the  /f2-norm  of  a  transfer  function  G(s)  =  C{sl  —  A)~^  B  subject  to 
CMC^  =  1.  The  optimal  C*  lies  in  the  direction  of  the  left  eigenvector  of  LcM~^  corresponding 
to  the  largest  (smallest)  eigenvalue  of  LcM~^  for  maximum  (minimum)  HG'H?,  where  Lc  is  the 
controllability  Gramian,  defined  by  ALc  +  LcA^  +  BB^  =  0. 

Remark  5.3:  For  optimal  disturbance  rejection,  when  the  pair  (C,  A)  in  (5.1)  is  not  observ¬ 
able,  Theorem  5.1  yields  a  choice  of  B-i  corresponding  to  a  null  eigenvector  of  Lo.  This  choice 
leads  to  the  invariance  of  the  output  with  respect  to  the  disturbance,  i.e.,  [CB^  C AB^ . . .  CA'^'^Bi] 
=  0. 

Remark  5.4:  For  unstable  SISO  systems,  one  can  introduce  the  notion  of  a  a-shifted 
norm  defined  as  follows: 


IIGIlL  =  dt  , 

—  CLcaC^  =  B^LoaB  ,  Q  >  0, 

where  a  is  chosen  so  that  (A  —  al)  is  Hurwitz,  and  Lea  and  Loa  are  the  controllability  and 
observability  Gramians  of  the  pair  (A  —  oel^B)  and  (C,  A  —  al)  respectively.  Here,  the  optimal 
B*  would  then  lie  along  the  direction  of  the  eigenvector  of  M~^  Loa  corresponding  to  the  largest 
(smallest)  eigenvalue  of  M~^Loa  for  maximum  (minimum)  ||G|l2o- 

5.5.3  QU  ALITATIVE  PROPERTIES  OF  OPEN  LOOP  SYSTEMS  WITH  THE 
//2-OPTIMAL  ZEROS 

This  section  examines  several  control-theoretic  properties  of  systems  with  the  input  vector 
chosen  according  to  Theorem  5.1.  The  dual  results  for  the  optimal  measurement  vector  C  can 
be  obtained  similarly. 

5. 5.3.1  Optimal  Control  Effort 

The  optimal  transmission  from  u  to  y  is  formalized  in  Section  5.5.2  as  the  /f2-optimization 
of  the  corresponding  transfer  function.  A  question  arises  :  What  does  £f2-optimization  of  G(s) 
means  in  terms  of  the  control  effort  necessary  to  accomplish  a  given  task.  An  answer  is  given 
below. 

Consider  the  system, 


i(t)  =  Aa:(<)  -b  Bu{t)  ,  (5.5) 

y(0  =  Cz\t) .  xenr,n,yeiR 

Let  U  be  the  set  of  all  bounded  measurabl  controls  u{t)  that  transfer  (5.5)  from  the  initial 
condition  x(0)  =  Zo  to  a  final  condition  satis  ig  y{U)  =  y\  during  [0,ti].  Define 

fti 

J{xo,yutx)=  min  u^(T)dT  . 
u(t)eu  Jo 

Theorem  5.2:  Assume  that  A  in  (5.5)  is  Hurwitz  and  M  in  (5.2)  is  positive  definite  and 
symmetric.  Then  min  lim  J(xo,yi,fi)  {  max  lim  */(xo, yi, <i) )  is  achieved  at  R  which  is 
collinear  with  the  eigenvector  of  Lo  corresponding  to  its  largest  (smallest)  eigenvalue. 
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Proof:  First  we  prove  that 


J{^o^  y\^i\)  = 


_ (yi  _ 

/o' ' 


Indeed,  a  control  u°(t)  that  transfers  (5.5)  from  the  initial  condition  j(0)  =  to  a  final 
condition  satisfying  y(<i)  =  j/i  over  [0,ti|  can  be  calculated  as  follows; 


u\i)  =  -  Ce-^'-xo],  t  €  [0,<i]  , 


where 


JO 

This  control  is  optimal  in  the  sense  that  if  u^{t)  is  any  other  control  in  U,  then 

{u\t)f  dt>  [u'‘(t))'‘dt  . 


To  prove  this  we  note  that 


Therefore, 


Then 


yi  -  Ce^‘'xo  =  r  Ce^^''-''^Bu^{T)dT  =  /“  . 

Jo  Jo 

r  -  u°(T))dr  =  0  . 

Jo 

[yi  -  Ce^^xo]  -  u°(T))dr  =  0  , 


which  is  equivalent  to 


I*'  {u\t)  -  u°(t))  (0,^i)  [yi  -  Ce'“'xo]  dr  =  0  . 


From  the  last  expression  and  (5.7),  we  obtain 


Therefore, 


J  (^u^(r)  —  u^irfj  u°(r)dr  =  0  . 

J  (u'(r))  dr  =  J  ^u^(r)  —  u°(t)  +  u°(t)^  dr 

+2  («* (r)  -  u°(t))  u°{T)dT 

=  /  (u*(7-)-u°(r))^dr  +  ^‘  (u°(r))^dr 


which  results  in  (5.8). 
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Thus,  it  follows  from  (5.8)  that  the  minimum  energy  required  to  transfer  the  system  (5.5) 
from  i(0)  =  lo  to  a  final  condition  satisfying  r/(<i)  =  j/i  is  indeed  given  by  (5.6).  Hence  if  .4  is 
Hurwitz, 


lim  J{xo,yi,ti)  = 


Since  B^LoB,  subject  to  the  constraint  B^MB  =  1,  achieves  its  maximum  (minimum)  at  B 
collinear  with  the  maximum  (minimum)  eigenvalue  of  M~^ Lo,  Theorem  5.2  is  proved.  Q.E.D. 

Theorem  5.2  offers  an  alternative  characterization  of  /f2-optimal  transmission:  The  optimal 
Bi  (Bj)  ensures  minimum  control  (maximum  disturbance)  energy  in  transferring  the  system 
state  from  an  arbitrary  initial  state  to  a  final  state,  consistent  with  the  specified  output,  over  an 
asymptotically  infinite  time  interval. 

5. 5. 3. 2  Controllability  Properties 

As  stated  in  Theorem  5.1,  system  (5.1)  with  Bi=^Bi*  has  i/j-optimal  (maximal)  transmission 
properties  from  u  to  y.  A  question  arises:  Is  there  always  at  least  some  transmission  from  u  to 
every  component  of  x  in  such  a  system?  An  answer  is  given  below. 

Theorem  5.3:  Assume  A  in  (5.1)  is  Hurwitz,  M  in  (5.2)  is  positive  definite  and  symmetric 
and  the  maximum  eigenvalue  of  the  matrix  M~^Lo,  where  Lo  is  the  observability  Gramian  of 
(5.1),  is  simple.  Then  system  (5.1)  with  Bi=B\*  and  B2  =  0  is  controllable  if  and  only  if  there 
exists  no  left  eigenvector  of  A  orthogonal  to  the  eigenvector  of  M~^Lo  corresponding  to  its  largest 
eigenvalue. 

Proof:  Follows  directly  from  Theorem  5.1  and  the  PBH  controllability  test  [94].  Q.E.D. 
Theorem  5.4:  Assume  A  in  (5.1)  is  Hurwitz  and  B\=Bi*,  B2  =  0.  Assume  further  that 
M  =  I,  the  observability  Gramian  Lo  has  distinct  eigenvalues  and  rankr  1  =  n.  Then 

J 

system  (5.1)  is  controllable  only  if  {Lo,LoA^  —  A^Lo)  is  observable. 

Proof;  The  theorem  is  proved  by  contradiction.  Suppose  jBj  =  Bj  and  B2  =  0.  If 
(lo,  LoA^  -  A^Lo)  is  unobservable,  then  for  some  A,  there  exists  x  ^  0  such  that 

LoA^  -  A^Lo  J  ^  ’ 

i.e.  _  LoA^x  =  XA^x  ,  A  ^  Amax  • 

Since  A  is  a  simple  eigenvalue  of  Lo  with  a  unique  eigenvector  x,  this  implies  that  AJx  must 
be  collinear  with  x,  i.e.  there  must  exist  p  such  that  A^x  =  px.  This  shows  that  x^  is  a  left 
eigenvector  of  A.  Since  both,  x  and  B\*,  are  eigenvectors  of  Lo, 

x^b;  =  0, 

x^A  =  px^  . 

Therefore,  according  to  the  PBH  test,  system  (5.1)  with  By=B\*  and  B2  =  0  is  uncontrollable 
[95].  ^  Q.E.D. 

5. 5. 3. 3  Coordinate  Independence 

Theorem  5.1  may  yield  optimal  B\  and  B2  which  result  in  coordinate  dependent  zeros.  The 
following  statement  gives  condition  on  M  which  ensures  coordinate  independence. 

Theorem  5.5:  Assume  that  A  in  (5.1)  is  Hurwitz  and  M  in  (5.2)  is  positive  definite  and 

symmetric.  Let  5  be  a  similarity  transformation  and  M  be  the  representation  of  M  in  the  new 
basis.  Then  system  (5.1)  with  Bi=B*  and  52=5^  yields  zeros  which  are  independent  of  the 
coordinate  transformation  5  if  M  =  S^MS. 

Proof:  From  Theorem  5.1,  B*  and  Bj  satisfy 


M'^LoBl 

m-^LoB; 
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Under  the  similarity  transformation  5,  system  (5.1)  has  the  following  representation; 

A  =  S-^AS,  Bi  =  =  =  CS  and  L,  =  S^LoS  . 

The  optimal  B^  and  are  given  by  M~^LoBi  =  Amax^j  and  M~^  LoBl  =  respectively, 

i.e., 

M-^S'^LoSBl  =  AmaxA*  and  Ar'S^LoSB;  =  \^\^B;. 

Hence,  if  M  =  5^M5, 

b;  =  s-^Bi  and  b;  =  s-'^b;  . 

Since  optimal  B^  and  ^.re  related  to  B^  and  respectively,  by  the  transformation  5,  the 
zeros  defined  by  B^  and  are  independent  of  the  coordinate  transformation  5.  Q.E.D. 

If  M  =  C^C  +  A^C^CA  +  ...  4-  {A^-'^fC'^CA--\  then  M  =  5^(C^C  +  A^C^CA  + 
...  +  {A'^~^Y C A^~^)S  —  S^MS.  Therefore  this  normalization  indeed  results  in  coordinate 
independent  zeros.  If,  however,  M  =  I,  the  resulting  B2-optimal  zeros  may  be  coordinate 
dependent.  Indeed,  consider  the  system  (5.1)  with  M  =  I  and  A,C  as  follows: 


J  J]  .c  =  [ioi. 


In  this  case,  Bj  is  given  by 


b;  = 


so  that  the  zero  of  the  system  is  at  0.  Introducing  the  similarity  transformation 


results  in 


-3  -2 
1  0 


,  C  =  [0  1)  . 


This  gives 


and  the  zero  lies  at  -3.277.  Hence,  the  normalization  M  =  I  may  lead  to  coordinate  dependent 
zeros. 

5.5.3.4  Non-minimum  Phase  Properties 

It  should  be  pointed  out  that  the  transfer  functions  Gtt(s)  and  Grv{s)  associated  with  the 
optimal  vectors  BJ  and  Bj  may  not  be  minimum  phase. 

Theorem  5.6:  Assume  M  =  I  and  the  maximum  and  minimum  eigenvalues  of  the 
observability  Gramian  Lo  are  different.  Suppose  (5.1)  is  in  the  observer  canonical  form  with 
Bi  =  BJ  and  B2=BJ.  Then  at  least  one  of  the  transfer  functions,  Gu(s)  or  (7u,(s),  is  non¬ 
minimum  phase. 

Proof:  For  system  (5.1)  in  the  observer  canonical  form,  the  elements  of  vectors  Bj  and  Bj 
correspond  to  the  numerator  coefficients  of  the  transfer  functions  Gu(5)  and  G„;(s)  respectively. 
Thus,  Gtt(s)  and  G«,(s)  are  minimum  phase  only  if  all  the  elements  of  Bf  ana  Bj  are  of  the 
same  sign.  We  show  below  that  this  is,  in  fact,  impossible.  Indeed,  Bj  and  Bj  are  collinear  with 
two  eigenvectors  of  the  observability  Gramian  Lg.  As  a  symmetric  matrix,  Lo  cannot  have  two 
eigenvectors  with  all  components  of  the  same  sign.  To  prove  this,  assume  that  x  and  y  are  two 
eigenvectors  of  Lo<,  associated  with  two  eigenvalues,  A  and  fi  {fi  ^  \)  respectively.  Assume  that 
X  has  all  components  of  the  same  sign.  Then,  since  x  and  y  are  orthogonal,  i.e.  Xtyi  =  0, 
not  all  components  of  y  are  of  the  same  sign.  Thus,  Bj  and  BJ  cannot  have  all  components  of 
the  same  sign.  Q.E.D. 
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5.5.4  //2-OPTIMAL  zeros  IN  CLOSED  LOOP  ENVIRONMENT 


5. 5. 4.1  Problem  Formulation 

In  the  closed  loop  environment,  a  number  of  problems  concerning  ff2-optinial  zeros  can  be 
formulated.  The  following  three  are  of  importance  in  various  control-theoretic  situations; 

5. 5. 4. 1.1  /f2-optimal  zeros  in  the  disturbance  rejection  problem 

Here  the  system  (5.1)  is  considered  in  a  closed  loop  form,  for  instance,  with  u  =  Kx.  Define 
G{s)  =  C{sl  —  A  —  B\K)B2  and  we  want  to  choose  B^  such  that  IIGII2  is  minimized  subject  to 
IIB2II  =  1  (or  B2MB2  =  1,  A/  =  >  0).  Obviously,  if  (A  +  BiK)  is  Hurwitz,  the  answer  to 

this  question  follows  from  Section  5.5.2  above:  The  optimal  B2  is  collinear  with  the  eigenvector 
of  the  closed  loop  observability  Gramian,  Lf ,  corresponding  to  its  smallest  eigenvalue.  Here  Ll‘ 
is  defined  by 

{A  4-  BJ<fLl^  +  Ll^{A  +  BtK)  +  C^C  =  0. 

Due  to  the  separation  property,  the  same  result  holds  not  only  for  state  feedback,  but  also  for 
observer- based  output  feedback. 

5. 5. 4. 1.2  /f2-op^inial  zeros  in  the  regulator  problem 
Consider  the  system 


X  =  Ax  -h  Bu  , 

(5.10) 

y  =  Cx  ,  X  e  IR^,u,y  e  IR 

E{xo}  =  0,  E{xoxl}  =  Q  , 

where  Xo  —  z(0).  Assume  that 

u  =  Kx  , 

(5.11) 

where  K  is  chosen  so  that 

J  E  (x^CxCix  4-  dt|  , 

(5.12) 

is  minimized.  If  {A,B)  is  stabilizable,  (Ci,/!)  is  detectable,  it  is  well  known  that  the  optimal  K 
is  given  by 

K  =  -B'^P  ,  (5.13) 

where  P  is  the  positive  semi-definite  solution  of 

A^P-hPA-PRR^P-t-CfCi  =0  .  (5.14) 

The  optimad  value  attzdned  by  functional  (5.12)  with  feedback  (5.13)  is  : 

r  =  rr(PC?)  .  (5.15) 

Obviously,  the  value  J*  may  be  further  improved  by  minimizing  J*  over  all  admissible  vectors 
B.  To  cast  this  problem  in  H2-noTTn  minimization  form,  we  define  Q  =  BiPf  and  rewrite  (5.14) 

^  {A-BB^P/2fP  +  P{A-BB^P/2)  +  C^Ci=0.  (5.16) 

Hence 

r  =  Tr{PBiB()  =  Tt{B,PBI)  ,  (5.17) 

and  from  (5.16)  and  (5.17),  the  functional  J*  is  equivalent  to  the  square  of  the  /f2-norm  of  the 
following  transfer  function 


Gj{s)  =  CxisI  -A-  BKl2)-^Bx  . 


(5.18) 
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It  should  be  noted  that  the  matrix  (A  +  BI\f2)  is  Hurwitz  due  to  the  fact  that  LQ  designs  have 
a  guaranteed  gain  margin  from  1/2  to  -foo  [96].  Also,  for  a  given  A,  B\  and  Ci,  the  choice  of  B. 
which  determines  K,  will  generally  affect  both  the  zeros  and  poles  of  the  transfer  function  Gj{s) 
(see  Example  5.2). 

Thus,  the  problem  of  H2-optimal  zero  placement  regulator  problem  can  be  stated  as  follows: 
Given  (5.10),  (5.11),  (5.13)  and  (5.14),  find  vector  B  satisfying  the  constraint  ||B||  =  1  (or 
B^MB  =  1,  A/  =  M^,M  >  0)  such  that  is  minimized. 

A  necessary  condition  for  the  solution  of  this  problem  is  given  in  Section  5. 5. 4. 2. 

5.5.4. 1.3  i/2-optimal  zeros  in  the  servomechanism  problem 
Consider  again  system  (5.10)  and  assume  that 

u  =  Kx  +  V  ,  (5.19) 

where  v  is  an  exogenous  signal.  The  control  gain  K  is  assumed  to  be  chosen  according  to  the 
designer’s  favorite  methodology,  such  as  pole  placement,  LQG,  Hoo  design,  or  any  other.  The 
problem  of  /f2-optimal  zero  placement  in  the  servomechanism  problem  is  formulated  as  follows: 

Given  (5.10)  and  a  design  methodology  for  K,  find  B  satisfying  \\B\\  =  1  (or  B^MB  =  1, 
M  =  M  >  0)  which  maximizes  |1G„||2  where 

G^{s)  =  C{sI  -  A- BK)-^B  .  (5.20) 

A  solution  to  this  problem  is  given  below  under  the  assumption  that  K  is  designed  using  the 
LQR  approach,  i.e., 

K  =  , 

A^P  A  PA-  PBB'^P  +  GfCi  =  0  . 

Remark  5.5:  As  it  follows  from  (5.18)  and  (5.20),  the  main  difference  between  the  open  and 
closed  loop  ^f2-optimal  zero  placement  is  in  the  fact  that  in  the  latter  case,  vector  B  affects  not 
only  the  zeros  but  also  the  poles  of  the  trzmsfer  function. 

Remark  5.6:  The  problems  formulated  above  can  be  extended  to  dynamic  output  feedback 
as  well.  Indeed,  assume  that  we  use  a  controller 

It  =  Kx  +  V  ,  (5.21) 

X  =  Ax  A  Bu  +  L{]j  —  Cx)  .  (5.22) 

Then,  it  is  well  known  that  the  closed  loop  transfer  function  from  u  to  y  is  also  given  by  (5.20). 

5.5.4. 2  Main  Result 

5.5.4.2.1  /f2-optimal  zeros  in  the  regulator  problem 

Theorem  5.7:  Consider  the  closed  loop  system  (5.10),  (5.11),  (5.13)  and  (5.14).  Assume 
that  (i4,R,  Cl)  is  stabilizable  and  detectable.  Then  ^min  ||u’j||2  is  attained  at  vector  B  only  if 

it  is  collinear  with  an  eigenvector  of  matrix  S  =  PGP,  i.e., 

(PGP  -  \1)B  =  0  ,  (5.23) 


where  P,G  £  are  defined  by 

{A  -  BB^PfP  A  P{A  -  BB'^P)  +  PBB'^P  +  C/'C,  =  0  ,  (5.24) 

{A  -  BB'^P)G  +  G(A  -  BB'^Pf  +  PiPf  =  0  .  (5.25) 


IIGjII"  =  a  +  CiGCf  . 


In  addition, 
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Proof:  The  Lagrangian  for  the  problem  is 

L(B,P,G,^,X)  =  iBjPBx  +  Tt{G[A^P  +  PA  -  PBB^P  +  CfC,])  +  \{B'^ B  -  1), 

where  7,  A  €  IR  and  G  €  G  =  G^,  are  the  Lagrange  multipliers  which  are  not  simultanc 

ously  zero.  The  first  order  necessary  conditions  for  an  extremum  are  as  follows: 

=  -fBiBj  +  AG  +  GA^  -  GPBB^  -  BB'^PG  =  0  , 

=  -2PGPB  +  2AB  =  0  . 

Since  A  —  BB^P  is  stable  with  (Gi,A)  detectable,  7  =  0  implies  A  =  0,  G  =  0  in  the  above 
equations.  Hence,  it  can  be  assumed  without  loss  of  generality  that  7  =  1.  Furthermore,  G  is 
nonnegative  definite. 

Equations  (5.23)  and  (5.25)  are  thus  obtained  and  (5.24)  is  equivalent  to  (5.14).  From  (5.25), 
we  have 


dL 

dP 

dL 

dB 


Tr{PB,Bj) 


-Tr{P[A  -  BB'^P]G)  -  Tr[PG\A  -  BB'^P]^) 
-2Tr(P(A  -  BB^PJG) 

-2rr(BAG)  +  2Tr{PBB^PG) 

^2Tr{PAG)  +  2B^BGBB  . 


FVom  (5.23),  (5.24)  and  knowing  B^B  =  1,  we  have 

IIGjII^  =  TriPByBf) 

=  -2Tr{PAG)  +  2X 

—  — A  +  C\GCj  +  2A 
=  X  +  CjGCf  . 


Q.E.D. 

5.5.4.2.2  Bj-optimal  zeros  in  the  servomechanism  problem 

Theorem  5.8:  Consider  the  closed  loop  system  (5.10)^,  f5.19),  (5.13)  and  (5.14).  Assume 
that  (A,B,Gi)  is  stabilizable  and  detectable.  Then  max  IIGvIlz  is  attained  at  vector  B  only  if 

it  is  collinear  with  an  eigenvector  of  matrix  S  =  —  PGP  —  2PLfL^^  i.e.. 


(Z?  -  PGP  -  2PLiL^  -  XI)B  =  0  , 

where  P,L‘J,G  and  €  JR"*"  are  defined  by  the  following  equations: 

{A-BB^PfP  +  PiA-BB^P)  +  PBB^P  +  C^Ci  =  0, 
{A-BB'^PfL^  +  L^{A-BB^P)  +  C^C  =  0, 
(A  -  BB^B)G  +  G(A  -  BB^P)^  -  LfLl^BB^  -  BB^I^'l^'  =  0, 
(A  -  BB^P)Lf  A  LfiA-  BB^Pf  +  BB"^  =  0. 


(5.26) 


(5.27) 

(5.28) 

(5.29) 

(5.30) 


In  addition, 

max  IIGJl^  =  CLl‘C'^  . 
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Proof:  By  forming  the  Lagrangian. 

L{B,L^,P,G,Ll‘,‘f,X)  =  •,B'^Ll‘B-¥Tr{G[A'^P-^-PA-PBB'^P  +  CjC,]^ 

UjliA  -  BB'^PfL‘J  +  -  BB'^P)  +  C^C])  - 

where  7, A  €  iR  and  G  =  G^,Ll‘  =  (L^)^  are  the  Lagrange  multipliers  which  are  not  simul¬ 
taneously  zero.  The  closed  loop  observability  Gramian  Ll‘  is  defined  by  (5.28).  The  following 
conditions  are  necessary  for  an  extremum; 

^  =  AG  +  GA^  -  GPBB"^  -  BB'^PG  -  L‘JLI‘BB'^  -  BB'^LI‘L‘J  =  0  , 

(5.31) 


^  =  iA-BB^P)Ll‘  +  Ll‘iA-BB^Pf  +  yBB^  =  0, 

II  =  2'rL^B-2PGPB-4PL^L^B-2XB  =  0  . 


(5.32) 

(5.33) 


Since  A  —  BB^P  is  stable  with  {C\,A)  detectable,  7  =  0  implies  =  G  =  0  and  A  =  0  in 
(5.32),  (5  31)  and  (5.33)  respectively.  Hence,  it  can  be  assumed  without  loss  of  generality  that 
7  =  1.  Furthermore,  Zj  is  nonnegative  definite. 

Equations  (5.26),  (5.29)-(5.30)  are  thus  obtained  and  (5.27)  is  equivalent  to  (5.14).  From 
(5.28)  and  (5.30)  ,  we  nave 


TriBB'^L^)  =  -2Tr{AL'^L^)  +  2Tr{L^PBB^L'^) 
=  TriLfC'*'C)  =  CL^C^. 


Q.E.D. 

Remark  5.7:  Since  equations  (5.23)-(5.25)  and  (5.26)-(5.30)  are  all  coupled,  they  are 
formidable,  indeed.  Their  relative  utility  may  be  justifiea  by  the  following  two  arguments.  First, 
these  equations  can  be  used  to  verify  whether  vector  B  optimal  for  an  open  loop  system  remains 
optimal  r\  the  closed  loop  environment.  Section  5.5.5  below  examines  this  relationship.  Sec¬ 
ondly,  iterative  procedures  can  be  suggested  for  solving  these  equations  (see  below).  Although 
the  convergence  conditions  for  these  procedures  are  presently  unknown,  it  is  easy  to  show  that 
they  always  generate  sequences  which  contain  convergent  subsequences.  Indeed,  since  the  opti¬ 
mization  is  performed  over  the  compact  set  {B  :  B^B  =1},  any  sequence  contains  at  least  one 
convergent  subsequence  (97). 

The  following  iterative  procedure  was  tested  on  several  examples  for  solving  equations  (5.26)- 
(5.30): 


1)  Choose  a  vector  B{n),  for  n  =  0  such  that  B^{n)B{n)  =  1  (e.g.  B{n)  =  [100  ...  0]^). 

2)  Solve  the  algebraic  Ricatti  equation  (5.27)  for  n  =  0  giving  P(n). 


3)  Solve  the  Lyapunov  equations  (5.28)  and  (5.30)  for  n  =  0  giving  Lf{n)  and  L^‘{n). 

4)  Solve  the  Lyapunov  equation  (5.29)  for  n  =  0  giving  G(n). 
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5)  Solve  the  eigenvalue  equation  (5.26)  for  n  =  1  such  that  A(n)  is  the  largest  real  eigenvalue 
of  L‘J  —  PGP  —  2PLl‘Ll‘  and  B{n)  is  the  corresponding  eigenvector. 

6)  Iterate  steps  (2),  (3),  (4),  and  (5)  for  n  =  1,2,3,...,  giving  sequences  {P(n)),  {Ll‘{n)}, 
{i;'(n)),{G(n)),{A(n)},and 

A  similar  iterative  procedure  for  solving  equations  (5.23)-(5.25)  is  obtained  by  making  the 
appropriate  modifications  on  the  above  algorithm. 

5.5.5  RELATIONSHIP  BETWEEN  THE  OPEN  LOOP  AND  CLOSED  LOOP  H2- 
OPTIMAL  ZEROS 

The  /f2-optimal  input  vector  B  calculated  according  to  Theorem  5  .1  for  the  open  loop 
environment  may  or  may  not  be  optimal  in  the  closed  loop  ctise  pertaining  to  Theorem  5.7. 
Below  we  illustrate  this  assertion  by  two  examples  and  formulate  a  conjecture  as  to  when  one  or 
another  situation  may  take  place. 

Example  5.1:  Consider  the  system  (5.10)  with  A  and  C  as  follows: 


[0  -10000  ■ 

^  [  1  -200  ’ 

C  =  (Olj. 


In  this  case,  the  open  loop  observability  Gramian  is 


(5.34) 


Lo 


2.5E  -07  o' 

0  2.5B  -  03  J  ’ 


and,  therefore,  vector  B",  that  achieves  max  |(G||2,  where  (^(s)  =  C(s/  —  A)~^B,  is  given  by 

l|B||=i 


Next,  substituting  (5.34)  in  equations  (5.23)-(5.25),  we  calculate  the  vector  B  optimal  in  the 
regulator  problem.  Let  Ci  =  C,Bi  =  [1  1]^,  and  M  =  I.  Invoking  the  iterative  procedure 
described  in  Section  5.5. 4.2,  we  obtain 


■  25.0122  0  ]  [  2.5E  -  07  0  ‘ 

0  0.0025  J  ’  ^  [0  2.5E  -  03  J  ’ 


and,  finally. 


which  was  verified  by  performing  an  exhaustive  search  over  the  set  {B  :  ||B||  =  1}.  Thus,  in  this 
particular  situation,  the  optimal  open  loop  B*j  and  the  optimal  closed  loop  B^  are  the  same. 
Example  5.2:  Consider  again  (5.10)  with  A  and  C  given  by 


C  =  |0  1), 


lA 


and  Cl  =  C,M  =  I.  The  open  loop  optimal  for  this  case  turns  out  to  be 


whereais  the  regulator  problem  results  in 


^  rO.4384 
10.8988 J  ’ 

K  =  (0.0347  0.1465]  . 

Obviously,  the  two  vectors  5*/  and  differ  significantly.  Note  also  that  the  zero  of  the  triple 
{A,  Bi,  Cl)  is  at  -1  while  that  of  the  triple  (.4  +  B*iK/2,  Bj ,  Cj)  lies  at  -0.984. 

Based  on  our  experience  in  computing  the  optimal  input  vector  B  in  open  loop  (Theorem 
5.1)  and  closed  loop  (Theorem  5.7),  we  are  led  to  the  following  conjecture. 

Conjecture  5.1:  When  the  LQR  methodology  is  used  to  compute  the  feedback  gain  K  in 
the  closed  loop  case,  optimal  vectors  B*;  and  B“;  are  close  to  each  other  for  system  with  “fast” 
open-loop  poles. 

The  reason  behind  such  a  conjecture  is  that  for  systems  possessing  open  loop  poles  in  the  far 
left  half  complex  plane,  the  location  of  the  closed  loop  poles  will  not  be  significantly  different 
from  that  of  the  open  loop  poles  with  a  LQR  controller.  Hence  optimal  vector  B*;  obtained  from 
the  open  loop  system  dynamics  will  be  close  to  that  of  B*/. 

However,  in  general,  in  a  closed  loop  system  with  a  specified  feedback  co’-  ‘  -oiler,  for  maximum 
Bj-norm  transmission  from  an  exogeneous  input  to  the  output  y,  the  choici.  of  the  optimal  input 
vector  must  take  into  <  onsideration  the  dynamics  of  the  closed- loop  system  and  not  the  open-loop 
dynamics.  To  be  more  precise,  the  choice  of  the  input  vector  B  and  the  controller  design  should 
be  performed  simultaneously  as  presented  in  Theorems  5.7  and  5.8. 
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9  APPENDICES 


APPENDIX  Al:  LEMMAS  A.1-A.5 


Lemma  A.l:  Let 


PW  =  E  A, A,  , 


where  A  =  {A, . A^}  is  a  set  of  self  conjugate  points  in  D{a.u).  The  function  P  achieves  its 

global  majcimum  subject  to  the  constraint 

m 

A.  =  a  , 

is  I 

at  Xi  ^  dDyi  =  1, . . .  ,m,  where  dD  denotes  the  boundary  of  D{(t,u). 

Proof:  Consider  a  set  of  self  conjugate  points  A  =  {Aj, . . . ,  Am).  Then 


m  mJ2 

S  XiXj  =  5I(Re  Ai)^  +  (Imag  A,)^  +  4  ^  (Re  A,)(Re  A>)  , 

•j«l  isl 

•o  •<) 


(A.l) 


where 


jlmag  A,j  <  —Re  A,  . 


Prom  (A.l)  and  (A.2),  it  follows  that  maxA6D(<r*;)  P{^)  >s  achieved  at 


|Imag  A,1  =  —Re  A,-  , 


i  —  1, . . . ,  m  f 


(A.2) 


i.e.  when  A*  €  dD,  i  —  1, . . . ,  m.  Q.E.D. 

Consider  the  optimization  problem:  For  m  even,  and  cr^  —  u;*  >  0,  find  a  self  conjugate  set  of 
eigenvalues  A  =  {Ai, . . . ,  A„j}  €  D{<t,u)  such  that 


P(x) = x;  . 


(A.3) 


is  maximized  subject  to 


X)  A,  =  a  <  0 

is] 


(A.4) 


Lemma  A.2:  The  global  solution  A*  =  {AJ, . . . ,  AJ^}  of  problem  (A.3),  (A.4)  is  given  by 


A‘  =  — (<T+jw)  , 
ma 


(A.5) 

(A.6) 
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t  =  l,...,m/2  , 
i  =  l,...,m/2  , 


and 


(A.7) 


Proof:  As  a  consequence  of  Lemma  A.l,  the  global  solution  A  of  (A. 3).  subject  to  the 
constraint  (A.4),  must  have  the  form 


A.  =  r,(<r  +  ja;)  ,  r,  >  1  ,  2  =:l,...,m/2 

'^m/2+.  =  ri(<r-ju;)  ,  j  =  l . m/2 


Therefore, 


Tft 

p(A)  = 

tjml 

•<J 

‘•<J 

and  the  constraint  (A.4)  reduces  to 


m/2 


Since  <r^  —  >  0,maxP(A)  is  achieved  at  r  which  maximizes  J2T/mi  r,Tj. 


•<j 


(A.8) 


(A.9) 


Define  the  set  r  =  {rj, . . .  Since  the  constraint  (A.9)  is  always  regular,  we  apply  the 

Lagrange  multiplier  rule.  The  Lagrangrian  is 


m/2 

/m/7 

^  riVj  - 

1,^.1 

•<} 

\'=1 

conditions 

dL{rJ) 

=  0, 

dr 

dL{r,0) 

=  0, 

d0 

a 

2cr 


yield  a  linear  system  of  equations  whose  solution  is 


r,  a: . . .  _  — 


mcr 


together  with 


2- 
2<7m 


The  second  order  conditions  ensure  that  (A.5)-(A.6)  is  a  strict  maximum.  Mo-eover.  since  under 
the  constraint  (A.9)  the  cost  function  (A, 8)  is 
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quadratic,  this  maximum  is  the  global  maximum,  and  the  proof  is  complete. 

Q.E.D. 

Consider  the  optimization  problem  :  For  m  even,  and  <  0,  find  a  self  conjugate  set 

of  eigenvalues  A  =  {A] , . . . ,  A„}  €  D{a,j;)  such  that 

P(\)  =  f;  A, A,  ,  (A.IO) 

•.^-1 

»<} 

is  maximized  subject  to 

f;A,=a.  (A.U) 

tsl 


by 


Lemma  A.3;  The  global  solution  A*  =  {Aj,...,A;j}  of  problem  (A.IO),  (A. 11)  is  given 


K  =  r-ia  +  ju). 


1  =  1,.. .  ,m/2  , 
I  “  1 , . . . ,  mjO,  , 


where 
In  addition. 


a  fm  \ 

'■i  =  ^  -  (-2  “  1  j  ^  *  =  2 . m/2 


(A.12) 

(A.13) 

(A.I4) 


P(A-)  =  (I  -  l)  (f  -  i)  +  •  (A.15) 


Proof;  As  a  consequence  of  Lemma  A.l,  the  global  solution  A  of  (A.IO),  subject  to  the 
constraint  (A.  11),  must  have  the  form 

A,  =  r,(o  +  >a;)  ,  r,  >  1  ,  t  =  l,...,m/2, 

Am/2+<  =  ri{<T-ju))  ,  t  =  1, . . . ,  m/2  . 


•<j 


Since  o’  —  w*  <  0,  and  from  (A. 8),  max  P(A)  is  achieved  at  r  which  minimizes  53" Hence, 
we  need  to  determine  the  set  r  =  {rj  •  <  •  r^/z}  which  minimizes 

m/2 

^(’■)  =  51  riVj  , 


•o 


subject  to  the  constraints 


r.  >  1  , 

m/2 

^  a  m 

ir  2<''2 


t  —  If,..,  Tn/2  f 


(A.16) 

(A.17) 


Since  the  constraints  (A, 16),  (A.17)  are  regular,  we  apply  the  Karush-Kuhn-Tucker  rule.  The 
Lagrangian  is 


m/2 


m/2 


/m/2 


L{r,  fij)=  52  ~  H  (»•.  -  1)  -  I  H  ^  ^  - 


•  .i-l 
•<; 


1=1 


i»=l 


86 


where  the  non-negative  Lagrange  multiplier  fi  »  satisfies  /i,(r,  —  1)  =  0.  for  all 

*  =  The  first  order  necessary  conditions  yield  a  linear  system  of  equations  whose 

solution  generates  a  set  of  candidate  extremal  points  given  by 


~  —  k 

r-  =  2z _ 1  u 

2 

r.  =  1  ,  fii 

0  = 

2c  =  -  t  ’ 


A*.  =  1  , 

//,  =  0  , 


i  ^  v4(r)  . 
«  €  /l(r)  , 


(A.IS) 


where  i4(r)  =  {r  :  r  =  1}  is  the  set  of  active  constraints  and  k  <  m  denotes  the  number  of  active 
constraints.  It  can  be  shown  that 

Hence,  among  the  candidate  extremal  points  given  by  (A. 18),  the  solution  k  =  m/2  -  1  yields 
the  global  minimum 


together  with 


.  a 

»•;  =  1, 


^  =  f-. 

..  =  0,  .  =  ^-1 


t  =  2,..., m/2  , 


t  =5  2, . . . ,  m/2  . 


The  results  (A.12)-(A.15)  then  follow  immediately. 


Q.E.D. 


Consider  the  optimization  problem:  For  m  even,  and  >  0,  find  a  set  of  self  conjugate 

eigenvalues  A  =  {Ai, . . . ,  Am}  €  D{a,u})  such  that 


P(A)  =  •£  A,Ai  , 


(A.19) 


is  maximized  subject  to 


Aq  =  ToO  , 

m 

^  A,  =  a  . 


ro>l  , 


(A.20) 

(A.21) 


Lemma  A.4:  The  global  solution  A*  =  {Ai,...,AJ^}  of  problem  (A.19)-(A.2i)  is  given 


..  a  -  ro(r  . 

Ai  =  — — + 
m<T 

\  •  O  “  I  ^ 

=  ——  (<^  -  3^)  . 
mtr 
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i  =  l,...,m/2  , 
i  =  l,...,m/2  , 


and 

P(A”)  =  roo-(a  -  rocr)  +  ((m  -  1)<T^ • 

Proof:  From  (A. 20),  the  expression  (A. 19)  can  be  expressed  as: 

P(X)  =  , 

•o 

m 

=  ro<7(fl  -  raff)  +  ^  A, Ay  . 

(A.22) 

•<» 

Also  the  constraints  (A. 20)  and  (A.21)  are  equivalent  to  the 
constrzunt: 

following  single 

m 

^  A,  =  a  -  roa  . 

i=l 

(A.23) 

Since  the  first  term  in  the  right  hand  side  of  (A.22)  is  constant,  the  solution  of  the  optimization 
problem  (A.22),  (A.23)  is  obtained  from  Lemma  A.2.  Q.E.D. 

Consider  the  optimization  problem  :  For  m  even,  and  <t^  <  0,  find  a  set  of  self  conjugate 

eigenvalues  A  =  {Ai, . . . ,  A,„}  €  Di<T,u))  such  that 

P(A)  =  £  A.A,  , 

(A.24) 

is  maximized  subject  to 

Ao  =  roo  ,  ro  ^  1  , 

(A.25) 

II 

(A.26) 

Lemma  A. 5:  The  global  solution  A*  =  {A^, . . . ,  AJ^}  is  given  by 

A-  =  r,.(<7  +  ju>),  i=l,...,m/2, 

Am/2+,  =  ^■  (<^  -  jw)  ,  X  =  1, . . . ,  m/2  , 

(A.27) 

(A.28) 

where 

rl=  (2  1)  “d  ’’i  =  1  >  1  =  2,..., m/2 

(A.29) 

In  addition. 

P(A  )  =  ro<T(a  roo-)  +  (  2  )  ( 2  er  ) 

(A.30) 

Proof:  The  results  (A.27)-(A,30)  follow  immediately  from  (A. 22),  (A. 23),  and  Lemma 
A.3.  Q.E.D. 
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APPENDIX  A2:  DERIVATION  OF  AVERAGED  EQUATIONS  FOR 

HELICOPTER  WITH  VHHC 


A2.1  AVERAGED  EQUATIONS  FOR  HELICOPTER  IN  HC'.  ER 

Let  ii  =  i,i2  =  ^1^:3  =  x,X4  =  6.  Equations  (4.21)  and  (4.22)  are  equivalent  to 


■  Xl 

X3 

.  *4  . 

- 1 _ _  -  *■> 

/(m|  —ffi)  sin^  ) 


-mj/rJ  CMZ2  - /kxi -f  ao*inwl  -  (X3  I 

—mjigunxj  +  ptinNii/t  -  tfZi  J 


=  jr,(».x)+i.Vj  (i.x)  . 

The  generating  equation  [20]  for  this  case  is  : 


)**"(f) 


«(mi^m3*in9  «))  V</ 

_ raiiSJU-^— «in 


olinT 

(mi  +m2~mj#in*X2 ) 
~tt«inx^  tinr 
qini+mj— m»»«f?X2) 


where  t  sst/c. 

We  obtain  the  general  solution  of  the  above  equation  in  the  fast  time  r  as 


—a  cos  r 

,  .  2  ’ 

mi  +  m2  —  m2  sin  C2 

a  sin  X2  cos  r 

- T - - - r-2— +  C4  , 

mi  +  m2  —  m2  sin  C2 


where  c,,  t  =  1, . . . ,  4,  are  constants. 

Conser  .ntly,  the  substitution  for  this  case  is 


Z2 

*»OrCO«T 


mi  «in^  xj 

g  in  cot  r 


I-  ^  J  L  <(mi+mj— mjMn^  *3)  ^  J 

and  the  equation  in  Bogoliubov’s  standard  form  C2ui  be  written  as  [20] 


■ISl' 

f» 

t  t  0 


Xxltghir.z)) 


2«m>>  >m  »<»  c»t  r 

*a)^ 

^2®22aLS!!Lii2,£2t*l£*li  ^ 

(mi  ^ mf  ein^  *3) 


^  iiii  *3) 

ft  ft  _^___-****  ^***^  — _.— 

^  I(m|^m9«-m3  •in^  wil/*("*l+***a*""**  **•  •*) 

-a  €••_*- _ ^ 

m|  m3  mh^  *3 

_ a  ofi  ca*  r  , 


Ma<*  Jj) 

(^^.;r,:c»;».,) + *0*  “•**  ■  *'*  ■"  ‘  ^ 

-ma<«un23  +0»mNui  -  n  7,)  +  *♦) 


m3»in^  <3 


Therefore,  applying  the  averaging  principle  ,  we  obtain  equations  (4,23)  and  (4.24)  where 
IT  =  z  and  15  = 


A2.2  AVERAGED  EQUATIONS  FOR  HELICOPTER  IN  FORWARD  FLIGHT 

Let  a:,  =  x,X2  =  y,X3  =  e,XA  =  i,xs  =  y,*6  =  Equations  (4.4)-(4.6)  are  equivalent  to  ; 


1  0  0 
0  1  0 
0  0  1 
0  0  0 
0  0  0 
0  0  0 

*4 

*•  + 
-ma&tj  co« *3  -  tel  +  ao  <»•»*•«»>  wJ  -  (»4 
m3<^  unxs  +  (mi  m,)p  unit  -  tei  Ubm  +  <*0  imuMaw*  -  *»*» 
-mtttmnmt  +  finnNwt-iist 


The  generating  equation  for  this  case  is  : 


where  t  ss  t/t. 

The  general  solution  of  the  above  equation  is  : 

*1  =  Cl  , 

X2  =  C2  , 

X3  =  C3  , 
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>1 


X4  - 

Xs  = 

X6  - 


-a[ 
“( 


(mi  +  m;  —  m;  cos*  C3)  cos  /i  +  m^  sin  C3  cos  C3  sin  fi 

mi(mi  +  m2) 

(mi  +  m2  —  m2  sin^  C3)  sin  ^  +  m2  sin  C3  cos  C3  cos  n 


sin  C3  cos  fi 


mi(mi  +  m2) 
cos  C3  sin  fi ' 


mi£ 


cos  r  -I-  Ce  . 


where  c,-,  z  =  1, . . . ,  6,  are  constants. 
The  resulting  substitution  is 


■  Xx  ■ 

X2 

X3 

X4 

Xi 

.  ^6  . 

Zl 

Z2 

23 


_  /  (mi  -f mt-ma  co«^  «3 )  eoa  n-f  «»n *»  co» »a  tin  u  \  T  +  C 
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_  /  (mi-fm}-m}  »in^  »3)«iim4-m?  liuxs  co>tj  cot(i\  T  +  2 

“  \  mi(mi+m2)  /  '"  ® 

/  «in  n  < 


and  the  equation  in  Bogoliubov’s  standard  form  is  [20] 


dt 

dr 


s  c 


where 


r  1 
0 
0 

0 

0 

0 


1  0 
0  1 
0  0 
0  0 
0  0 
0  0 
0 


0 

0 

1 


(c-  «.  cm  M-«i«  »»  «» I.)  ^  ^ 


mj 


0 

0 

0 

1 

0 

0 


I 

0 

0 

0 

0 


0 

0 

0 


m|(mi  +m3) 
•!»  €94  tn 

fa 


mi* 


0 

0 

0 


0 
0 
0 

mT^ 

—  co»  t\ 
mjt 

mTmaT*  J 


!(mi  »a>€o»M^ma  »»■  •a  co»  la  tin 

(mi 4m7— "«»«■*  iiCO»»aco«  m 

mi(mi+m») 

a  ("«  co»r  +  H 

+aoco.^MDu;l -Ol 

m,/a2  ^)’«n*3  +  (m,  +  m,)*  UAM  -  fc*j  Uirn  +  ao«n4«nu.t  - 

-mjtg  .in <3  +  <3.inNw»  -  n  (!■*»  +  „) 


((mi  +  m2  —  m2  cos^  23)  cos  /i  +  m2  sin  23  cos  23  sin  ^l)a  cos  t 
'  ~  mi(mi  +  m2) 

((mi  +  m2  —  m2  sin^  23)  sin  fx  +  m^  sin  23  cos  23  cos  n)a  cos  t 


<1  =  za  - 
<2  =  25  — 


mi(mi  +  m2) 


Upon  applying  the  averaging  principle,  we  ''Main  equations  (4.31)  where  21  =  2,22 
zS  =  <t>. 


V  and 


91 


